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Hi Derivers, 

The Texas Instruments website has a free, collaborative forum where you can now share 
your favorite Derive activities and upload Derive worksheets and math files. The online 
Activities Exchange allows educators to post an activity that has proven effective in class 
- or browse submissions from colleagues. 

For more information on the Activities Exchange, visit 


http://education.ti.com/exchanqe . 

From this page, you can browse by subject or TI product; search for specific activities; or 
submit a new activity. When you click on the "Submit an activity" link you will be asked 
to log onto the TI web site. (If you don't already have a User Name and Password, sim- 
ply register to create one.) Afterwards, you will be guided through a series of screens to 
input the activity information and upload activity file(s). Please complete all the online 
template components (e.g., "before the activity, "during the activity", "other materials", 
"after the activity", etc.) as more information encourages activity usage. 

Below are links to example Derive activities: 

http://education.ti.com/educationportal/activityexchange/activity_detail.do7activityid =41 
70&cid=us 

http://education.ti.com/educationportal/activityexchange/activity_detail.do7activityid =41 
72&cid=us 

Regards, 

Theresa Shelby 
Texas Instruments 


I am very grateful for permission to present one sample lesson in this DNL as a mouthwatering (page 
33), Josef 


CAME 2005 Symposium: Second Announcement 
Shaping Research and Development of 
Computer Algebra in Mathematics Education 

To be held in October 19-20, 2005 in Roanoke, Virginia USA 
in conjunction with PME-NA-27, October 20-23, 2005 
http://www.pmena.org/2005/ 

Themes of the Symposium 

1. CAS , instrumentation and the anthropological approach 
Presenter. John Monaghan, University of Leeds, UK 

2. The impact of CAS on our understanding of mathematics education 
Presenter. Werner Peschek, University of Klagenfurt, Austria 

3. Teacher learning while teaching with CAS 

Presenter. Baerbel Barzel, Universitat Duisburg-Essen, Germany 

Each of the three themes will be addressed in plenary lectures and in topic groups, where the issues 
will be worked on in more detail, based on the experiences of the participants. The topic groups will 
aim at producing a synthesis of the discussions in the plenary and topic group sessions. In the discus- 
sion groups there will be, apart from time set aside to discuss issues arising from the plenary ad- 
dresses, opportunity for further short presentations. 

The plenary papers will appear on CAME’s Website before the Symposium: 

http : /_/www . lonklab . ac . uk/came/ 
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Dear DUG members, 

When I started collecting stuff for DNL#58 I 
intended to end up with about 40-44 pages. As you 
can see it became more than 50 pages. There is one 
contribution from Rudeger Baumann not only an- 
swering an until now unanswered challenge but 
also offering a fine introduction in treating strings 
and in programming with them. Speaking about 
programming Titbits 30 must be mentioned. Johann 
Wiesenbauer recommends his Titbits because of 
containing some examples of very efficient pro- 
gramming. In my opinion Titbits 30 dealing with 
elliptic curves are one of his very best in his long 
series of Titbits published in the Newsletters so far, 
many thanks Johann, outdoing yourself once more. 

I am very grateful for permission to reprint 
Carl Leinbach's paper on GPS which can be found 
among the many activities which can be 
downloaded from the Tl-website (Information 
page). CarTs contribution is a mouthwatering ex- 
ample standing for many other papers. There is 
only little experience needed to transfer the activi- 
ties from PC to handheld and vice versa. 

Please notice the many announcements of use- 
ful websites and conferences presented in this 
DNL. 


Last weekend I was in the happy position to 
attend US AC AS 3 in Atlanta. There is a wonderful 
group of enthusiastic US-teachers led by Jim 
Schultz, Natalie Jakucyn and Bob Me Collum who 
want to propagate the use of CAS in US curricu- 
lum. 

Natalie Jakucyn made it possible to organize a 
CAS-strand at the 2006 NCTM Regional Confer- 
ence in Chicago. Please support this great idea by 
submitting a CAS-related paper (deadline 11 July 
2005, see last page! !). 

Thanks also to Ernest Carpenter from South 
Africa - he urged to publish the Coons article - for 
his mails and the provided websites on “NURB- 
surfaces”. 

I wish you a wonderful summer and am looking 
forward to providing DNL#59 and the revised 
DNL#7 in fall. 

Best regards to you and your families 



Some Australian websites: (University of Melbourne) 

extranet . edf ac . unimelb . edu . au/DSME/CAS-CAT/ 
extranet . edf ac . unimelb . edu . au/DSME/RITEMATHS/ 
edf ac . unimelb . edu . au/DSME/ 

www . vcaa .vic.edu.au/ vce/ studies/mathematics/caspilolt/ cas index . html 

2003 PME Proceedings online: online . terc . edu/ 

Canadian website: www . edu . gov . on . ca 

Individuals on Ontario schools who could provide additional information: 

Thomas_Steinke6occdsb . on . ca and f red. f erneyhough6peelsb . com. 

A Derive dedicated website in Belgium: http : / /cage . ugent . be/~svw/DfW/ 


Download all DNL-DERIVE- and Tl-flles from 


http : / / www . austromath . at /dug/ 
http : // www .bk- teachware . com/main . asp?session=375059 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & CAS-77 User Group. 
It is published at least four times a year 
with a contents of 44 pages minimum. The 
goals of the DNL are to enable the ex- 
change of experiences made with DERIVE 
and the 77-89 /92/ Titanium! Voyage 200 as 
well as to create a group to discuss the 
possibilities of new methodical and didac- 
tical manners in teaching mathematics. 

As many of the DERIVE Users are also 
using the CAS-77v the DNL tries to com- 
bine the applications of these modem tech- 
nologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone/FAX: 43-(0)2275/8207 
e-mail: nojo.boehm@pgv.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & CAS-77 Newsletter will be. 


Next issue: September 2005 

Deadline 15 August 2005 


Preview: Contributions waiting to be published (selection) 

Pringles, B. Grabinger, GER 

Two Stage Least Squares, M. R. Phillips, USA 

Some simulations of Random Experiments, J. Bohm, AUT & L. Kopp, GER 

Wonderful World of Pedal Curves, J. Bohm 

Another Task for End Examination, J. Lechner, AUT 

Tools for 3D-Problems, P. Luke-Rosendahl, GER 

ANOVA with DERIVE & Tl, M. R. Phillips, USA 

Hill-Encription, J. Bohm 

Farey Sequences on the Tl, M. Lesmes-Acosta, COL 
Simulating a Graphing Calculator in DERIVE , J. Bohm, AUT 
Henon & Co, J. Bohm 

Challenges from Fermat, Bj. Felsager, DEN 
Are all Bodies falling equally fast, J. Lechner, AUT 
Modelling Traffic Density, Th. Himmelbauer, AUT 

Do you know this? Cabri & CAS on PC and Handheld, W. Wegscheider, AUT 
An Interesting Problem with a Triangle, P. Luke-Rosendahl, GER 
Diophantine Polynomials, Duncan McDougall, CAN 
Rosettes, J. Lechner, AUT 
Jacobi-lteration, Heinz Rainer Geyer, GER 

and Setif, FRA; Vermeylen, BEL; Leinbach, USA; Koller, AUT; Baumann, GER; 
Keunecke, GER and others 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wurmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 
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Steve Schonefeld 


sschonef eld6mchsi . com 


Hi Joseph & Noor. 

Every time I look at my poster from the Montreal DERIVE/TI conference, I remember the folks in 
charge of that conference confiding that the graphics for that poster were not created using DERIVE. 
They were created using MAPLE. I hope we can remedy that for the next conference. DERIVE can 
create many beautiful and interesting graphics. Perhaps members of DUG can contribute graphics for 
the next poster. Let me start the ball rolling with the attached contribution. It is a torus wrapped in a 
torus knot. 

Wishing you and yours a happy Easter. 

Steve Schonefeld 



Great idea! 

So let's open a contest for the DES-TIME 2006 
poster. 



DES-TIME 2006 in Dresden, Juli 20 - 23, 2006 


www.des-time-2006.de 
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Tania Koller, St.Polten. Austria 

Some of Tania's students had problems with annotating new plots in the 2D Plot Window. 

I came across this strange behaviour, too and I could not find the reason for its occurrence. 

I wanted to plot y = x, y = x A 2/4+x/2-1 and f(x):=2 sin(3x) with annotating the graphs, but the variables 
didn't appear as shown on the screen shot. 

Any idea what happened?? 



Theresa Shelby 

Thank you very much for the information and the file. The mystery is solved! 

After the annotation is created, Derive substitutes the Horizontal Axis Title for every occurrence of x 
and the Vertical Axis Title for every occurrence of y in the annotation of a plotted expression. Cur- 
rently, if a title in the plot window's Options > Display > Axes is blank, then an empty string will be 
substituted into the annotation. So, in this case, if both the horizontal and vertical axis titles are blank, 
"y=x" becomes "=" and "y = x A 2/4+x/2-1" becomes "= A 2/4+/2-1". 

In order to keep the horizontal and vertical titles invisible, yet still preserving the annotations, I advise 
using the Options > Display > Axes tab to restore the default titles of "x" and "y" and choose a font 
color for the axes titles that's close enough to the background of the plot window (e.g., "white") so as 
not to be visible. 

For the next release of Derive, I will modify the code so that it does no substitution of the Axes Titles in 
the annotation if the title is empty. 


Aloha, 

Theresa 
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Ernest Carpenter (South Africa) 

Hi Josef, 

Thanks for the files relating to the Coons-Surfaces. I enjoyed it tremendously (especially after it was 
translated©) for mainly the following reasons: 

> Firstly it was clear and nicely presented and made a topic, that usually make me feel more 
math-illiterate than I’m currently is due to the sometimes complex representation in literature, 
a little bit more understanding! 

> Working with CAD (Microstation, AutoCAD, etc.) and modeling systems (Amapi, Maya, 3D 
Max 5, Cinema 4D, etc.), I constantly use these NURB related surfaces without always fully 
understanding the main principles behind the construction - now it makes more sense! 

> It prompted me to ‘research’ the tools I use and get a better understanding of how big and 
complex curve and surface modeling are! (Interesting history - 

http://rt001 kxn.eresmas.net/principal/history CAGD.pdf) 


(one of the packages I also use is Eovia’s Amapi and a nice tutorial showing how some of these sur- 
faces are used in modeling is: http://server2.eovia.com/tutorials/amapi6/shoe/) 


Right is a quick CAD representation of a NURB surface 
using a degenerated Coons-Gordon algorithm (Microsta- 
tion) 

I also immensely enjoyed the new and revised DNL’s. 

As an exercise I ‘modeled’ the “Two Underground Tunnels” 
problem in CAD to see how the geometric results compare 
with the numerical results in DNL05. They are the same©! 
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Rudeger submitted this contribution some time ago, but this does not influence its inter- 
esting nature. His article isiwritten in german and I will try to provide short English sum- 
maries. In the original paper all functions, variables and programs have German names. 

To make it easier for our non-German speaking members I use English names. There 
are two DERIVE-files which can be downloaded: engel.dfw with the German names and 
engel_eng.dfw with the English names. Most of the problems and solutions are self ex- 
planatory. I added the respective solutions for our many TI-CAS users. Josef 

Engelfolgen und Schieberegister 

Zum Umgang mit Zeichenketten und Zahlen 

Engel Sequences and Shift Registers 
How to treat Strings and Numbers 

Riideger Baumann, Celle 

Wie viele Anfragen und insbesondere auch die Leserzuschriften in DNL 49 zeigen, besteht nach wie 
vor Bedarf an elementaren Einfiihrungsaufgaben und Ubungen zum Programmieren in Derive. Zwar 
hat Josefs Buch hier Abhilfe geschaffen - aber andererseits sind natiirlich noch Fragen offen geblie- 
ben. Im folgenden sollen einige elementare Uberlegungen zur Behandlung von Zeichenketten und 
Zahlen angestellt werden, die erfahrungsgemaB dem Anfanger Schwierigkeiten bereitet. Als Anwen- 
dung werden zwei interessante Aufgaben gelost. 

There are many requests on introductory examples and exercises in programming in Derive. In this 
contribution some elementary considerations on the treatise of strings and numbers are presented. 
Then using the acquired knowledge two interesting problems will be solved. 


Elementare Operationen 

Zeichenketten (engl.: strings), also (aus Buchstaben, Satzzeichen, Ziffem usw. gebildete) Worter und 
Satze werden zwischen obenstehende ,,Gansefu6chen“ (engl.: double quotes) eingeschlossen. Die 
Operationen fur Zeichenketten und Listen (in Derive: Vektoren) sind eng miteinander verwandt: man 
kann die Listenoperationen auch auf Worter, Zahlen und Texte anwenden. Die drei Beispiele behan- 
deln die wechselseitige Umwandlung zwischen Zahlen, Zeichenketten und Listen. 


Beispiel 1: Eine gegebene Zahl soil in eine Zeichenkette verwandelt werden - und umgekehrt. 
Example 1: A given number shall be transformed into a string - and vice versa. 

Die Aufgabe Zahl in Wort ( NumberToWord ) lasst sich sehr einfach mit Hilfe der Funktion STRING 
losen: 


#1: NumberToWo rd(n) := STRING (n) 

#2: NumberToWo nd (3 . 14159) = [/, 314159, 100000] 


314159/100000] 


#3 : NumberToWo rd (-3 . 14159) 


314159 ] 


100000 


This is the way, the result 
appears in Exact Simplification. 


# 4 : 


NumberToWo nd (3 . 14159) 
Approximate #4 and #6 
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#5: 3.14159 


#6: NurnberToWo rd (-3 .14159) 

#7: [-, 3.14159] 

#8: NAME_T0_C0DES(-3. 14159) = [45, 51, 49, 52, 49, 53, 57, 47, 49, 48, 48, 48, 48, 48] 

#9: NurnberToWo rd2 (n) := APPROX (CODE S_TO_NAHE(NAHE_TO_CODES(n))) 

#10: NurnberToWo nd 2 (-3 .14159) = -3.14159 


Die Umkehrung Wort in Zahl ist nicht ganz so einfach, denn ein Wort kann Zeichen enthalten derart, 
dass es sich nicht als Zahl interpretieren lasst. Wir wandeln das Zahlwort zunachst - mit Hilfe der 
Funktion NAME_TO_CODES - in eine Liste von ASCII-Nummem um: 

The reverse operation WordToNumber is not so easy, because there might be characters which can- 
not be interpreted as a number. At first we transform the number-word into a list of ASCII-Codes. 

#11: Wo rdToASCIILi st(w) := NANIE_T0_CODES(w) 

#12: Wo ndToASCIILf st(-3140) = [45, 51, 49, 52, 48] 

#13: Wo rdToASCIILi st(3 .14) = [49, 53, 55, 47, 53, 48] 

Auf diese Liste wenden wir die Funktion CODES TO NAME an, die aus einer Liste von ASCII- 
Nummem eine Zahl macht. Sie weist allerdings Besonderheiten beziiglich der Null auf: 

C0DES_T0_NAME ( [ 48 ] ) = 0 
C0DES_T0_NAME([48, 51]) = 3 

Das heifit, die vordere Null wurde weggelassen. Aus diesem Verhalten entstehen zuweilen Komplika- 
tionen. Insgesamt erhalten wir 

#14: Wo rdTo N urn be r (w) := C0DES_T0_NAHE(NAHE_T0_C0DES(w)) 

#15: Wo rdTo Number (0123) = 123 

#16: WordToNumber (-1.23) = -1.23 

Das Argument ist ein String, daher unter " zu schreiben. 


As you can see the leading zero is omitted. This might cause sometimes problems. The argument is a 
string and must be written under It works with negative integers, but not with negative decimal num- 
bers. 


#17: WordToNumber (-314) + WondToNumber(314) = □ 

#18: WordToNumber (-3 ,14) + WordToNumber(3 .14) = -3.14 + 3.14 

The CAS-Tls provide two functions stri ng (number) and expr(string) for these purposes: 


f t i gibr afca fcbtherlpr gn I ofc 1 ean Up! 1 

■ stringt.nj ntw(.nj 


□one 

■ ntw(3. 14159} 


"3. 14159" 

■ ntwC -3. 14159} 


" "3. 14159" 

■ ^2^ ] 


"3*ji^2" 

■ expKw} + wtn(w} 


Done 

■ wtnC "0123"} 


123 

■ wtn( 11 -1.23 11 } 


-1.23 

wtn<" “1 _ 23 11 >1 

MAIN RAD AUTD 


FUNC 7/30 
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Beispiel 2: Einer natiirlichen Zahl soil die Liste ihrer Ziffem zugeordnet werden - und umgekehrt. 
Example 2: A natural number shall be assigned to the list of its digits - and vice versa. 

Bei einer Zahl n kann man die i-te Ziffer (von links) durch n SUB i (bzw. nii ) herausgreifen: 

NumberToList(n) VECTQR(n , i , DIM(n)) 

#19 : i 

#20: NumberToLi st(3141590) = [3, 1, 4, 1, 5, 9, 0] 

#21: NumberToLi st2(n) VECTOR (k - 48, k, NAM E_T0_C0 DE S ( n ) ) 

#22: NumberToLi st2 (3141590) = [3, 1, 4, 1, S, 9, 0] 

Even entering NumberToList(03l4l590) the leading zero does not appear in the Algebra Window. 

Die Umkehrung Liste in Zahl geht wie folgt: Wir definieren eine Funktion asc durch 

#23: asc(z) := FIRST ( NAM E_T0_C0DES(z)) 

#24: asc (3) = SI 

Die Funktion NAME_T0_C0DES liefert eine Fiste; ihr entnehmen wir mittels FIRST das erste (und einzi- 
ge) Element. 

#2S : ListToNumber(v) := C0DES_T0_NAME (VECTOR (FIRST (NAME_T0_C0DES(k)) , k, v)) 

#26: ListToNumber([0, 3, 1, 4, 0]) = 3140 

Es gibt auch eine rein arithmetische Losung: 

NumberToLi stl(n) 

If n < 0 
#27: [] 

ADDQIN(M0D(n , 10), NumberToLi stl(FL00R(n , 10))) 

#28: NumberToLi stl(3141S90) = [0, 9, S, 1, 4, 1, 3] 

Die Ziffem werden in umgekehrter Reihenfolge aufgelistet! Dies kann manchmal niitzlich sein. Wiin- 
schen wir uns die Ziffern in normaler Reihenfolge, hangen wir sie mittels APPEND hinten an. Da diese 
Funktion Listen verkettet, machen wir aus der Ziffer M0D(n , 10) durch EinschlieBen in eckige Klam- 
mern eine (einelementige) Liste. Die Funktion lautet: 

NumberToLi st2(n) := 

If n < 0 
#29: [] 

APPEND(NumberToLi st2(FL00R(n , 10)), [M0D(n, 10)]) 

#30: NumberToLi st2 (3141590) = [3, 1, 4, 1, S, 9, 0] 

NumberToLi stl returns the digits in reverse order which can be useful sometimes. To obtain the regu- 
lar order we use append instead of addoin. Take care for the brackets enclosing the digit MOD(n,lO) in 
NumberToLi st2 making this digit to a list containing only one element. This is necessary, because 
append works with lists only. 


f » Ifi 1 gibr alca fclother fprgn I ole 1 

ean Up! 1 

■ ntol(31452} 

£3 1 

4 5 2> 

d i 1 i ) , 




■ S H 

k= 1 ■ 

k . 10 diKlii- 

k J ■* ltonCli) 




Done 

■ lton(£3 

1 4 5 2IO 


31452 

■ 1 tonC £3 

1 0 G 2 

0» 

310020 

■ lton(£0 

3 1G0 

2 or) 

31002G 

lton<<0. 

3, 1,0,0, 2 

,0>>l 


MAIN 

RAD AUTD 

FUNC £1/30 


Wi— |i = LV^LI 

Y i i v Y i'M- Y F5 Y ^ :i-r 

sic |OthM'|PrgnI0|CL:ix:i^ 





Done 

■ ntol(31452) £3 1 4 

5 2y 

d i 1 i ) 

■ Z 1: 

k= 1 - 

i k 10 diri ^ li ^ - k ) lton(li) 

Done 

■ lton( £3 

1 4 5 2IO 

31452 

... ,dim<strinsr<num>}}+ntol<num} 

MfilN 

F;fi[i fiUTD FUNC 1£/£1 
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Beispiel 3: Einer gegebenen Zeichenkette soil die Liste ihrer Zeichen zugeordnet werden - und umge- 
kehrt. 

Example 3: A given string shall be assigned to the list of its charcters - and vice versa. 

Die Aufgabe Wort in Liste lost sich analog zu Beispiel 2: 


WbndToLi'st(w) := VECTOR (w , i, DIMM) 

#31: i 

#32: Wo ndToLi st(DERIVE & CAS-TI) = [D, E, R, I, V, E, ,4, , C, A, S, T, I] 

Fiir die Umkehrung Liste in Wort orientieren wir uns ebenfalls an Beispiel 2: 

Li stToWordO(v) will resolve the problem with the omitted leading zero(s). 

#33: ListToWond(v) := CODES_TO_NAME(APPEND(NAME_TO_CODES(v))) 

#34: Li stToWondC[D, N, L, #, 5, 8]) = DNL#58 


it i -toy fit y f::t y fh^ y fe y fst y "i 
|t IfllgebralCalclotherlPrgEiIolciean Up| | 

■lton(£3 1082 

0» 

310020 

■lton(tO 3100 

2 0>) 

310020 

■ seqCnicKword ? k ? 1), k. 

1 ? din(wor-d)) + uto ► 
Done 

■ wtol( "DERIDE & CfiS-TI 

£ II 0 II II £ II ll p« II ll J II 

11 ) 

ll || ll ll ^ II 

II II 

■ ltow({"D" "H" "L" 

"#" i.^.i 

"8"}) 

"DHLttSS" 

„,w<{ M D M - M N M - M L M - 1 



MAIN RAD AUTD 

FUNC EH/30 



fFiT^OY FT? 

| t j - — | Contra 1 
: ltow< 1 i 
: Func 

: Local i ? w 




'YF::tYfHtY FF Y Fb T 

|lYQ|Oar|Find...|Mode 


:For i,l,din(li} 
: w& 1 i [ l ] -*w 
: EndFor 
: w 

: EndFuncI 


] 


MAIN RAD AUTD FUNC 


Die Engelfolge / The Engel Sequence 

In der Zeitschrift Praxis der Mathematik stellte Arthur Engel im Jahr 1967 folgende Aufgabe: 
1967 Arthur Engel posed a problem in Praxis der Mathematik: 


„Es werden die zehn Ziffem z u z 2 , ..., zi 0 beliebig angenommen. Diese Ziffemfolge wird nach 
dem Gesetz z n = z n _i 0 + z n _ 9 (mod 10) mit n = 1 1, 12, 13, ... fortgesetzt. Man zeige, dass die Folge 
periodisch wird, ja sogar sofortperiodisch. Was kann man iiber die Periodenlange aussagen?“ 

Let z 1? z 2 , ..., z 10 ten arbitrary digits. This sequence is continued according to the following rule: 
z n = z n _ 10 + z n _ 9 (mod 1 0) with n = 1 1 , 1 2, 1 3. Show that this sequence is periodical and that there 
is no preperiod appearing. 

Beispiel/ Example : 1234567890 3579135793826048... 

Wir wollen hier - als Anwendung der Uberlegungen zu Zeichenketten und Zahlen - die Erzeugung 
der Folge programmieren. Zunachst hat man sich zu entscheiden, ob die Folge als Liste, Zahl oder 
Zeichenkette geschrieben werden soil. Es handelt sich um Ziffem, die modulo 10 addiert werden: dies 
spricht fiir die Verwendung von Zahlen. Andererseits miissten wir dann auf weit zuriickliegende Zif- 
fem einer (groBen) Zahl zugreifen, die Zahl um eine Stelle verlangem und die neue Ziffem als Einer- 
ziffer anfugen: dies scheinen schwierige arithmetische Operationen zu sein. Somit empfiehlt sich die 
Verwendung von Zeichenketten. 

Instead of working with numbers we recommend the use of strings because we have to refer to num- 
bers lying far back and form new numbers. This seems to cause difficult arithmetic operations. 
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Wir stellen zunachst die Lange k des gegebenen Worts (Ziffernfolge) fest und bauen dann eine Schlei- 
fe (LOOP) auf, die zur Bestimmung der i-ten Ziffer (i > k) um k bzw. um k - 1 Schritte zuriickgreift: 
Folge SUB (i - k) bzw. Folge SUB (i - k + 1). Da es sich um Zeichen handelt, miissen sie in 
Zahlen umgewandelt werden (siehe oben), die sich dann modulo 10 addieren lassen und rechts ange- 
hangt werden. Hier wirkt die Funktion APPEND auf Zeichenketten (von Derive als Listen aufgefasst). 

Um die aus der Folge herausgegriffenen Ziffem als Zahlen behandeln zu konnen, miissen wir die Auf- 
gabe Ziffer in Zahl losen: 

As we have to treat the digits as numbers we have to solve the problem DigitToNumber. 


#37: Digi tToNumber(z) := FIRST (NAM E_TG_C0DES(z)) - 4S 


#3S : Digi tToNun(iber(3) = 3 


Zuerst wird der Ziffer (hier: 3) ihre ASCII-Nummer (51) zugeordnet und sodann die der Null (48) 
abgezogen: 51 - 48 = 3. 


The program arthur(start , length) with first results: 


Arthur(start := "1234B67S90 11 , length := 100, seq, zl, z2) := 

Pnog 

seq := start 
k := DIM (start) 
i := k 
Loop 

#39: i :+ 1 

If i > length 
RETURN seq 

zl := Digi tToNumber(seq 4 .(i - k)) 

z 2 := Digi tToNumber(seq.j.(i - k +1)) 

seq := APPEND(seq , STRING (M0D(zl + z2 , 10))) 

#40: ArthurQ 

#41 : 12 34 56 7B9 03 5 7913 5 79 3B 26 04B 26 210B64 2 OBB 31S4 06 2B6149 246B 0475316 04B412B4 7642253 0213 064 7B 3 2 34 36 01515 5779 

#42: Arthur(B62, 20) = B624B624B624B6248624 

#43: Arthur(555 , 20) = 55500505550050555005 

#44: Arthur(26B , 50) = 26SS4620B2S00S0SBS66420626SS4620S2S00S0SSB66420626 

#45: Arthur(220, 50) = 220424660262SS06B6440S4S220424660262SS06S6440S4S22 

#46: Arthur (5000, 50) = 50005005505055550005005505055550005005505055550005 

Can you find the periods? 

Arthur(862 , 20) = 8624 8624862486248624 
Arthur(555 , 20) = 5550050 5550050555005 

Arthur(268 , 50) = 268846208280080888664206 26884620828008088866420626 
Arthur(220 , 50) = 220424660262880686440848 22042466026288068644084822 
Arthur(5000 , 50) = 500050055050555 50005005505055550005005505055550005 


Wir finden bei dreistelligen Anfangswortem die Periodenlangen 4 (1), 7 (1), 24 (5), 28 (1), 168 (5); in 
Klammern steht jeweils die Anzahl der Zyklen. Probe: l-4+l-7 + 5*24+l*28 + 5-168 = 999. 
Die Periodizitat mittels Programm zu erkennen, scheint schwierig. Leider ist das in DNL#47, S. 38 
gestellte Problem der Erkennung von Perioden noch nicht gelost. 
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Rudeger notes that the problem to detect periods has not been solved until now (posed in DNL#47). 
This was reason enough for me to make a try: 


periods? (num , per, k, i, j) : = 

Prog 

kul 

Loop 

If k = DIM (num) 

RETURN "no period detected!" 

per := num 4 ,[l k] 

#47: i' := k + 1 

Loop 

j := IF(M0D(i, k) = 0, k, M0D(i , k)) 

If num.|.i ± per.|.j exit 

i :+ 1 

If i > DIM (num) 

RETURN per 
k :+ 1 

#48 : pen ode? (Arthur (5000, 50)) = 500050055050555 

#49: DIM (periods? (Arthur(220, 50))) = 24 

This program works, but to keep it more general it would be fine to also find out if there is a preperio- 
dical part existing. So look at the following: 


period? ( num , per, k, f, j, v, num_) : = 

Prog 

num_ := num 
v := 0 
Loop 
v :+ 1 
k := 1 
Loop 

If k = DIM (num) exit 

per := numj,[l k] 

#50: i := k + 1 

Loop 

j := IF(M0D(i , k) = 0, k, M0D(i , k)) 

If numj,i t perj,j exit 

i :+ 1 

If i > DIM (num) 

RETURN ["Preperi'od" , "Period"; num_j,[l, 
k :+ 1 

num := REST(num) 

If DIM (num) = 1 

RETURN "no period detected!" 


v - 1], per] 


#51: period? (Arthur (268 , 50)) 


P repe ri od P e ri od 

268846 2 08 28 008 0888664 2 06 


Some other examples to illustrate peri od?(st ring): 
5 

#54: 

173 


#55 : 0 . 028901734104046242774566473988439306358381502890173 

#56 : period? (028901734104046242774566473988439306358381502890173) 

Preperi'od Period 

0289 01734104 046 24 2 7 74 5664 7 39884 39 3 06 3 58 3815 _ 


# 57 : 
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#S8 : period? (00251B759tIL87B9trL8759Qlfi759(ILS759trL8759QlJB759) 

P repe n' od P e ri od 


#B9 : 


002 5 


1B7B90 


#60: pe n' od 7 (AB ABC ABC ABC ) = 


P repe nod Period 
AB ABC 


#61 : pe ri od 7 (ABC ABC DABC DAC ) 


P repe r i od P e ri od 

AB CABCDABCDA _ 


#62: 

#63: 


#64: 


#6B : 


peri'od?(ABCABCDABCDAE) = no period detected! 

pe ri od 7 (12 34 B 26SS46 2 OS 28 008 0SSS664 2 06 26SS46 2 OS 2S 00S 0SSS664 2 06 26 ) 
P repe ri od P e ri od 

_ 1234B 26SS4620S2S00S0SSS664206 . 

pe ri od 7 (12 34 B 26SS46 2 OS 2S 00S 0SSS664 2 06 26SS46 2 OS 2S 00S 0SSS664 2 06 2 7 ) 


#66: no period detected! 

DIM ((period? (Arthur (123, BOO))) ) = 16S 
#69 : 2,2 

DIM ((period? (Arthur (1234, 2000))) ) = 1B60 

#70: 2,2 


We can finally count all appearing periodics. #119 - #123 shows the end of this procedure. Rudeger 
started with "001" giving 999 three digit starting numbers. I start with "000" resulting in 1000 Engel- 
Sequences with a three digit initialisation number. 


#119: p0 := VECTOR (DIM (pen ode? (Arthur (v, 200))), v, VECTOR (STRING (k) , k, 700, 999)) 
#120: total := APPEND(pl, p2 , p3 , p4 , pB , p6 , p7 , pS , p9) 


#121: 

[[1, 

4, 

7, 

24, 

28. 

168], 

VECTG R ( DIM ( SE L ECT (x = v, v, total)), x, [1, 4, 7, 24, 28, 168])] 


" 1 

4 

7 

24 

28 

168 " 


#122: 

. 1 

4 

7 

120 

28 

840 . 



#123: I([l, 4,7, 120, 28, 840]) = 1000 


|Ti7»iriT 1 

1^ Icont-rol | 

1/0 

r fh t 'i 

Uar 

fFind...TnodeT I 


it 17 m' fit y f::t y fht y fe y fst y "i 
|t I fllgebralCalclotherlPrgnlolciean Up| 

: arthurtstart, le) 

: Func 

: Loca 1 se 7 k , i ? z_ 1 ? z_2 
: start- -*se: dim<:start)-*k 
: For i ? k+1 ? le 
: exprtnidtse, i-k ? l))-*z_l 
: exprtnidtse, i-k+l ? l))*z_2 
: seS:stringtnodtz_l+z_2 7 10))*se 
: EndFor 
: se 

: EndFunc 

■ arthurC "362" , 20) 

11 86248624862486248024 11 

■ arthurC "555" , 20) 

11 55500505550050555005 11 

■ art-hurC 11 5000" ? 50) 

"500050055050555500050055050555500050(^ 

art hur< "GOOD 11 , 50> 

MAIN RAD AUTD 

FUNC BATT 

MAIN RAD AUTD 

FUNC 3/30 


I'll leave the program for period detecting with the Tl as a challenge for the reader. 
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Die Erzeugung der Folge konnen wir uns mittels eines sogenannten Schieberegisters vorstellen (siehe 
Engel 1991, S. 105 ff). Ein Schieberegister (engl.: shift register) besteht aus einer Aneinanderreihung 
von Speicherzellen, die so geschaltet sind, dass - auf einen Taktimpuls hin - der Inhalt jeder Zelle an 
die linke Nachbarzelle weitergereicht wird. Dazu gehort eine Ruckkopplungsschaltung , deren Eingan- 
ge mit bestimmten Zellen verbunden sind, und deren Ausgang in den Eingang des Schieberegisters 
miindet. Die zugehorige (Pseudo-) Zufallsfolge wird von den Zahlen gebildet, die aus der Zelle am 
linken Rand heraustreten. 



Bildl /Figure 1 
Shift Register 


We can imagine creating the sequence using a so called 
shift register (Engel, 1991). This consists of a sequence of 
cells, which are connected in such a way that at each clock 
pulse the content of each cell is transmitted to its left neigh- 
bour. A reaction coupling is connected with certain cells as 
input and gives an output which feeds the sequence from 
the right side. The respective sequence of (pseudo) random 
numbers is formed by the numbers which leave the register 
at its left end. 


Die einfachste Version ist ein linear riickgekoppeltes Schieberegister (engl.: linear feedback shift regi- 
ster): hier besteht die Ruckkopplungsschaltung aus einer XOR-Verknupfung (Exklusiv-Oder) der an- 
gezapften Zellen, was mathematisch der Addition modulo 2 entspricht. 

Die Folge der angezapften Zellen, das Anzapfmuster (engl.: tap sequence), stellen wir durch einen 
Null-Eins-Vektor dar (in Bild 1 ware dies [1, 1, 0, 0, 0, 0]). Dann kann das neue Bit mittels skalarer 
Multiplikation (in Bild 1: MOD([l, 0, 1, 1, 1, 0] • [1, 1, 0, 0, 0, 0], 2) = 1) bestimmt werden. 

The simplest version is a linear feedback shift register (LFSR) whose input is the XOR-operation (ex- 
clusive or) of some of its outputs. These output cells are called taps. The tap sequence is represented 
by a binary vector. Here it is [1 ,1 ,0,0, 0,0] because the first two cells are connected by XOR giving the 
next cell which is here 1 XOR 0 = 1. This new bit can be calculated by a dotproduct modulo 2 of the 
register content vector and the tap sequence vector MOD([1 , 0, 1 , 1 , 1 , 0] ■ [1 , 1 , 0, 0, 0, 0], 2) = 1 . 

The following program for producing a (pseudo) random binary sequence uses working with strings. 
Let's assume that we have already produced the sequence 10111011001 and we want to generate 
bit #1 2 (i = 1 2). We need a "window" [6, 7, 8, 9, 10, 1 1 ] which selects the last k = 6 bits. 

Das folgende Programm macht von einer wichtigen Moglichkeit der Bearbeitung von Zeichenketten 
Gebrauch: man kann aus einer Zeichenkette durch Vorgabe einer Liste mittels der Operation SUB ge- 
zielt Elemente herausgreifen. Es sei beispielsweise die Folge 10111011001 bereits gewonnen, und nun 
soil das zwolfte Bit erzeugt werden (i = 12). Dazu bilden wir das „Fenster“ [6, 7, 8, 9, 10, 1 1], mit 
dessen Hilfe die k = 6 letzten Bits herausgegriffen werden: 


#74: 

sequence := [1, 0, 1, 

i, i, 


i, 

1, 

#7S : 

[k := 6, i := 12] 





#76: 

window := VECTOR (W - 

k + j , 

: , 

o, 

k - 1) 

#77: 

sequence = [0, 

wf ndow 

1, 1, 

o, 

o, 

1] 


Werden sie mit dem Anzapfmuster skalar multipliziert, ergibt sich das neue, in den Eingang zuriickge- 
fiihrte Bit: 
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#78: pattern := [1, 1, 0, 0, 0, 0] 

M0D(pattern- sequence , 2) = 1 

#79 : wi ndow 


Damit bekommen wir das folgende Programm: 

#80: Di'gi'tToNuniber(z) := FIRST(NAME_TQ_CQDES(z)) - 48 

WQndToLi'st2(w) := VECT0R(Di'gi tTo Number (w ), i, DIM(w)) 

#81: i 

Li'stToWondOCv , w) 

Prog 
w := 1111 
Loop 

#82: If v = [] 

RETURN w 

w := APPEND(w, FIRST(v)) 
v := REST(v) 

sr(start := "1 0111 0" , pattern := 11 110000" , length := 69, seq , pat, window, i , k) := 
Prog 

seq := Wo ndToLi'st2 (start) 
pat := Wo rdToLi s t2 (pattern) 
k := DIM (start) 

-i := k 
#83: Loop 

i :+ 1 

If i > length 

RETURN Li stToWo rdO(seq) 
window := VECTOR (i - k + j, j, 0, k - 1) 
seq := APPEND(seq, [MOD (pat- seqj,wi ndow , 2)]) 

#84: sr() 

#SB : lammiOCniQlQl Ol mil OOOOQlOOOQnOOQlQlO Ol 1 1 1 0100011100100101101110 

We investigate the periodic behaviour: 

#86: peri'od7(sr()) 


#87: 


Preperiod Period 

1011101100110101 01 11111 0000010000110001010 01 1 1 1 0100011100100101 


DIM((peri'od7(sr())) ) = 63 

# 88 : 2,2 


Das Schieberegister liefert die maximale Periodenlange 2 6 - 1 = 63. Die Erzeugung von Zufallszahlen 
mittels Schieberegister spielt in der Kryptologie eine wichtige Rolle (siehe Schneier 1991, S. 373 ff). 

Another LFSR with period length 127: 

#91: DIM (pen ode7(sr (1110001, 1010101, 200))) = 127 

Note: 

It might make sense to let students produce a number of elements of a LFSR-generated sequence by 
hands first to understand the process - and then check it using the program, Josef 


For example: start with 1 ,0,1 ,0,1 ,1 and use the tap sequence 1 ,1 ,0,0,1 ,1 
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1 

0 

1 

0 

1 

1 

1 

1 


1 


1 

0 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

◄ 1 







1 0 


t 

0 

t 

0 





1 

1 



1 

1 

◄ 1 







0 1 


1 

0 

1 

0 

1 

1 

1 

1 

1 

1 

1 

1 

1 

0 

◄ 1 








0 1 


Compare with the Derive output (you don't see the quotes in the Algebra Window!): 

sr("ioioii", "noon", 12) = "loionnnoi" 

Let's first investigate its periodic behaviour and then perform the calculation on the CAS-TI: 


pen'ode?(sr (101011 , nOOU, 100)) 

im mnmmi mnnm nnnmminmmmnm n m 1 1 1 nnnn m i mil nm i nnmi 
DIM (peri ode? (srClOCLGCLl, nOGCLl, 100))) = 63 


it i -toy fit y f::t y fht y fe y fst y "i 
|t I fllgebralCalclothenPrgroIolciean Up| I 


it i -toy F£t 1 
|t Icont-rol 

r f3t'i 

I/O 

'fht'i 

Gar 

fFi F nd...lModeT 1 

■ sr< 11 1G1 1 1G" , 11 1 10GGG 11 , 5G) 
"1G111G11GG11G1G1G11111 100000100001 1G(^ 

■ srC 11 1Q1Q11 11 , 11 1 1QQ1 1 11 , 30) 

"101011111101 1Q1QQQ1QQQG1G1 1QG1 11 

■ srC" 10101111" , "11001100" , 20) 

"10101111111 1GGG1GG1G" 

: sr-tstart, pat, le) 

: Func 

: Local k, i , win, se 
: dipi(star-t)-*k 
: start-^-se 
:For i,k+l,le 
smicKse, i-k, k)+win 

: se&str-ingtmodtdotP^winKwin) , winKpat) ) 
, 2) )->se 
: EndFor 
: se 

: EndFunc 

sr<"101 01111" - "llOOllOO" .20> 

MAIN RAD AUTU 

FUNC 3/30 

MAIN RAD AUTD 

FUNC 


pl8 


Riideger Baumann: Engelfolgen & Schieberegister 


D-N-L#58 


Die Schornfolge / The Schorn Sequence 

In DN# 49 (S. 39) wurde Engels Aufgabe von Richard Schorn wie folgt abgeandert: 

In DNL#49 (p 39) Richard Schorn varied Engel's problem as follows: 

Eine Folge von Zahlen soli dadurch erzeugt werden, dass Einer- und Zehnerziffer (modulo 10) 
addiert und vorne angehangt werden, die Einerziffer aber entfallt. 

Generate a sequence of numbers by adding the last two digits of an integer (modulo 1 0) and put- 
ting the result at the beginning of the number, omit the last digit. 

Beispiel/ Example: 3718, 9371, 8937, 0893, 2089, ... 

Bei Engel geht es somit um eine Folge von Ziffern, bei Schorn um eine Folge von Zahlen - ein be- 
trachtlicher Unterschied. Dennoch hangen beide Folgen eng miteinander zusammen. 

Engel Sequence is a sequence of digits, Schorn sequence is a sequence of numbers - quite a diffe- 
rence, but both kinds of sequences show a tight connection. 

fi 1 lZeros(word) is an auxiliary function to show the leading zeros (as in 0893). 

ffllZenos(wond , n) := 

Pnog 

Loop 

#111: If DIM (word) < n 

word := APPEND("G M , word) 

RETURN word 

Richard (start := "3718 11 , pi'ckout := [2, 1], length := 10, k, n0, sum_, rseq) := 

Pnog 

k := DIM (start) 

#112: n0 := Wo rdToN umber (start) 

sum_ := Z(M0D(FL00R(n , 10 A (j - 1)), 10), j, pi'ckout) 

rseq := ITERATES(FL00R(n , 10) + 10* (k - 1) -M0D(sum_, 10), n, n0, length) 
rseq := VECTOR (fi'llZenos (STRING (k_) , k) , k_, rseq) 

#113: Ri'chandQ = [3718, 9371, 8937, 0893, 2089, 7208, 8720, 2872, 9287, B92S , 0592] 

#114: DIM (Richard (3718, [2, 1], »)) - 1 = I860 

Another example: Instead of taking the last two digits we take the first two ones. Derive recognizes 
6798174312 as starting number of a repeating cycle. 

Richard (1234567890, [10, 9], ») 

[1234567890, 3123456789, 4312345678, 7431234567, 1743123456, 8174312345, 9817431234, 7981743123, SfrJrfJfcEHM . 3679817431, 
9367981743, 2936798174, 1293679817, 3129367981, 4312936798, 7431293679, 1743129367, 8174312936, 9817431293, 7981743129, 
6798174312] 


#117: 

DIM (Richard (5000, 

[2, 

1], 


- 1 = 15 

#118 : 

DIM (Richard (2000, 

[2, 

1], 


- 1 = 312 

#119 : 

DIM (Richard (1000, 

[2, 

1], 


- 1 = I860 


Richard bestatigt also die oben festgestellte Periodizitat der Folge mit 5000 als erstem Element. Die 
Gleichung 1 • 15 + 2 • 312 + 6 • 1560 = 9999 festigt die Vermutung, dass es einen Zyklus mit 15, zwei 
Zyklen mit 312 und sechs Zyklen mit 1560 Elementen gibt. Was sagt Arthur dazu? 

Richard confirms periodicity of the sequence initialized by 5000. Equation 1x15+2x312+6x1560=9999 
consolidates the conjecture that there is one cycle with 15, two cycles with 312 and six cycles with 
1560 elements. What is Arthur's opinion? 
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Arthur("2000" , 322) = 

"20002002202422664820202222444688046840424466802482620882860046040644608068864640004004 

404844228640404444888266082680848822604864240664620082080288206026628280008008808688446 

280808888666422064260686644208628480228240064060466402042246460006006606266882460606666 

22284402842026228840624686044648002802082280408448282000200220" 


DIM(periode?(Arthur("2000" , 322))) = 312 

Aufgaben / Tasks 

1 . Welcher Zusammenhang besteht zwischen den Engel- und den Schornfolgen? Kann man insbeson- 
dere aus der Periodizitat der Schornfolgen auf die der Engelfolgen schlieBen? 

Which connection is between Engel sequences and Schorn sequences? Is it possible to deduce 
periodicity of Engel sequences from periodicity of Schorn sequences? 

2. Untersuchen Sie die vom Schieberegister mit der Anfangsbelegung 10101110011 und dem An- 
zapfmuster 10100000000 erzeugte Folge. Hat sie maximale Periodenlange? 

Investigate the sequence generated by a shift register with initial sequence 10101110011 and tap 
sequence 10100000000. Is this sequence of maximum period length? 

3. In der Folge 1983113835952... ist jede Ziffer ab der funften die mod-10-Summe der vier vorange- 
henden Ziffern. Enthalt die Folge das Wort (a) 1234, (b) 3269, (c) 5198, (d) nochmals 1983? Expe- 
rimentieren Sie mit verschiedenen Anfangswerten. Wie hangt die Periodenlange vom Anfangswert 
ab? 

In the sequence 1983113835952... each digit starting with the fifth is the mod 10 sum of its four 
predecessors. Das this sequence contain the words (a) 1234, (b) 3269, (c) 3269, (d) again 1983? 
Experiment with various initial values. How does period legth depend on the intial value? 

4. Angenommen, ein Kryptanalytiker hat die Teilfolge 00100110 eines linearen Vier-Bit-Schiebe- 
registers abgefangen und mochte den weiteren Verlauf der Folge voraussagen. Zu diesem Zweck 
sucht er das Anzapfmuster zu ermitteln. Helfen Sie ihm! 

Assumed a cryptoanalyst intercepted a partial sequence 00100110 of a linear four-bit-shift register 
and wants to predict the following elements. For this purpose he tries to detect the tap sequence. 
Can you help him? 

Literatur / References 

Bohm, J.: Programmieren in Derive. Hagenberg: Kutzler-Verlag, 2001 

Engel, A.: Eine periodische Dezimalzahl (Aufgabe P 302). In: Praxis der Mathematik 9 (1967), S. 166 
Engel, A.: Mathematisches Experimentieren mit dem PC. Stuttgart: Klett-Verlag, 1991 
Schneier, B.: Applied Cryptography. New York: Wiley, 1996 
Welke, St.: Detection of Periods. In: DNL 49, S. 26-28 


Useful related websites 

www. elektronik-kompendium . de/ sites/ dig/ 0210211. htm 
home . ecn . ab . ca/~ j savard/ crypto/ co04 0801 . htm 

www . newwave instruments . com/ resources/ articles/m_sequence_linear_f eedback 
shift register lfsr.htm 
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Coons [1, 2] -Surfaces with Derive and CAS-TI 

Josef Bohm 


We start constructing a surface above the unit square in the first quadrant which is bounded by four 
connected straight lines. Two lines are lying in the xz-plane andyz-plane, the other two lines are lying 
in the parallel verical planes. 

We have four vertices: P 0 o(0,0,-l), Pi 0 (l,0,5), Pn(l,l,l) and P 0 i(0,l,3). 

We designate the line in the xz-plane as gxO (points on it have variable x- values and y = 0 = constant: 
gxO(P 0 o,Pio). Analoguous definitions for gxl(P 0 i,Pn), gOy(P 0 o,Poi) andgly(Pio,Pn). 


DERIVE can help producing a sketch of the situation. 


(Left matrix produces the vertex points and the right 
one connects these points to a closed polygon in 3D 
space.) 


' [ 0 , 0 , - 1 ] " 
[I, 0, s] 
[I, 1, I] 

. [0, 1, 3] . 


' 0 

0 

-1 ' 


— 1 

0 

5 


— 1 

1 

1 


0 

1 

UP 


. o 

0 

-1 . 



z 



It is easy to produce the four lines in parameter form. The parameters are running from 0 to 1. We 
have to take care that opposite lines show the same orientation. It is no problem to name the parame- 
ters asx andy! 


gx0(z00, zOl, zlO, zll, x) := [x, 0, zOO + x*(zl0 - z00)] 
gxl(z00, zOl, zlO, zll, x) := [x, 1, zOl + x*(zll - zOl)] 


gQy(z00, zOl, zlO, zll, y) := [0, y, zOO + y*(z01 - z00)] 
gly(zQ0, zOl, zlO, zll, y) := [1, y, zlO + y*(zll - zlO)] 


Via [gx0(-l,3,5,l) ,gxl(-l, 3,5,1) ,g0y(-l, 3,5,1) ,gly(-l, 3, 5,1)] and Insert > Plot 
in the 3D-Plot Window we obtain the representations of the segments with boundaries 0 and 1 for the 
parameter. Don't forget to activate the box Apply parameters to rest of plot list. 


We can imagine two points X 0 and X\ on two opposite lines 
gxO and gxl which are generated by the same parameter 
value x. We transfer these two points and the connecting 
segment to the graph. This segment lies parallel to the 
yz-plane and needs parameter y for any constant x. 

Running x and y between 0 and 1 we receive a set of points, 
which form a surface in 3D-space. We found its parameter 
representation - with parameters x and y (which are a and t 
for DERIVE) 


y 
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[xl := [x , 0, zOO + x-(zl0 - zOO) ] , x2 := [x , 1, zOG_ + x- ( zll - zOCL)]] 
z(z00, zlO, zll, zOl) := xl + y- (x2 - xl) 
z(-l, 5, 1, 3) 

[x, y, 2-x-(3 - 4-y) + 4-y - 1] 


Inspecting the parameter form of the surface we find that the third component is nothing else than 
explicit representation of the z- values ot the surface: 

z(x,y) = -8xy + 6x + 4y-l mit 0<x,y<l 


In DERIVE notation this reads as 


IF(0 <x<lA0<y<l, 2x(3 - 4y) + 4y - 1, ?) 

So let's wrap this into one function depending on the z-values of the four corner points: 

bi"lin_unitz(z00, zlO, zll, zOCL) := 

IfO<x<lAO<y<l 
(z(z00, zlO, zll, z01))r3 
? 

bf lfn_um - tz(-l, B, 1, 3) 




Rewriting the third component - z(x,y) - shows that it might be interpreted as a product of vectors and 
a matrix in two ways: 

x*(y*(z00 - zOl - zlO + zll) - zOO + zlO) + y*(z01 - zOO) + zOO 

z00*x*y - z01*x*y - zlO*x*y + zll*x*y - z00*x + zlO*x - z00*y + z01*y + zOO 

(1 - y)*((l - x) *z00 + x*zlO) + y * ( ( 1 - x)*z01 + x*zll) 


zOO 

zOl ' 


5k 

i 

r— 1 

i 

. zlO 

zll . 


. y . 


[1 - x, x]. 

[x*(y*(zOO - zOl - zlO + zll) - zOO + zlO) + y*(z01 - zOO) + zOO] 
[1 - y, y]. 


' zOO 

zlO " 


' 1 - X ' 

. zOl 

Zll . 


x 


|[x-(y-(z00 - zOl - zlO + zll) - zOO + zlO) + y-(zQl - zOO) + z00]| 
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i 

■ bi 1 in_u( "1,5, 1,3) x-(6 - f 

■x(6-Sy) + 4y- 1 ■+ zl(x 7 y) 

B f x -(6 - B ■ y) + 4- y - 1 , x > 0 and 
\mdef ? else 

| 1,1 ■» bil» 

[HBBE 

3 ■ y) + 4 ■ y - 1 
Done 
x ^ 1 and ^ 

... <6-8#y>+4#y-l,undef >+E2<x,y> 

lDDN:' RAD EXACT 3D H/30 



This leads to another - and more compact - definition of the surface: 

bil-in_un-it(zOO , zlO, zll, zOl) := 

If 0<x<lA0<y<l 

([1 - X, x]-[z00, zOl; ZlO, Zll].[l - y; y])U 

7 

b-il-in_unit(-l, 5, 1, 3) 

First generalization 


The base figure shall no longer be only a square but a rectangle with xO < x < xl and yO <y <yl.We 
change the vertices to 

Poo(0,0,l), P 5 o(5 ,0,6), P 53 (5,3,-l) and P 03 (0,3,2). 


Parameter representation of the segments must be 
changed for accepting the default parameter boundaries 
0 and 1 to run from xO to xl and from yO to y 1 . 

The sketch shows a top view for (x0,y0) = (2,1) and 
(xl,yl) = (5,6). Default values for function 
bilin_rect() are (0,0) and (1,1). 

bilin_rect ( zOO , zaO , zab, zOb, x0:=0, y0:=0, xl:=l, yl:=l, zxy):= 
prog ( 

zxy : = ( [l - (x - x0)/(xl - xO), (x - x0)/(xl - xO)]* 
[zOO, zOb; zaO , zab]* 

[1 - (y - yO) / (yl - yO) ; (y - yO) / (yl - yO) ] ) 4 , 1 , 
IF ( xO <x<xlAy0<y< yl , zxy, ?) 



bi"lin_rect(-l, 5, 1 , 3) = IF(x < 1a y < 1a x > 0a y > 0 , 2-x-(3 - 4-y) + 4-y - 1 , ?) 


This function seems to work properly and we try to produce two parallel surfaces. 



IF 


x<7^y<5^x>2/,y>2, 


15 


0 
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Additional exercise: 

Find the boundary lines and scale in such a way that the range for the parameter is from 0 to 1 for all 
segments, (see DERIVE file). 

[zrxOfl, 2, 6, -1, 0, 0, 5, 3), zrxb(l ? 2, 6, -1, 0, 0, 5, 3), 

zry0(l, 2, 6, -1, 0, 0, 5, 3), zrya(l, 2, 6, -1, 0, 0, 5, 3)] 



It might be useful for students to repeat the relationship between the parameter form and explicit form 
of the surface: 


[z_00 := [0, 0, 1], z_0b := [0, 3, 2], Z_ab := [5, 3, -1], Z_a0 := [5, 0, 6]] 
[g_ := z_a0 + t(z_O0 - z_a0) , h_ : = z_ab + t(z_0b - z_ab) ] 

9_ + s(h_ - g_) 

[5 - 5 ■ t , 3s, s ■ (8 ■ t - 7) - 5 ■ t + 6] 

S0LUE(x =5-5-t A y =3 -Sa z__ = s(8t - 7) - 5t + 6, [s, t, z__]) 

y 5 - x 5 ■ (y + 3) - x(8y - 15) 

s = ^ t = ^ z = 

3 5 15 


5- (y + 3) - x- ( 8 ■ y 


IF 


5 ■ Cy + 3) - x ■ (8 ■ y - 15) 

x<5 A y<3^x>0 A y>0, , ? 

15 


It is not really difficult to transfer this to the TI-89/TI-92/Voyage 200: 


rat 


. fT^ y f:=t y fht t f? y Y i 

|fllgebra|Calc|Qiher|Prgn)IQ|Clean Up| 1 


■ soloe(x = 5 - 5 ■ t- and y = 3 s and z = s(S^ 


^ ~(x ~ 5) 


and z = 


-(x (3 ■ y - 15) - 5 (y + 3)) 


5 “■ 15 

~(x(8y - 15) - 5(y + 5)) 


15 


^undef ? else 


? x > G and 

Done 


... >-5^(v+3))/15 J .undef )->zl(x J .v) 

MAIN RAD EXACT 3D E/30 
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Second generalization 

The second generalization is much more challenging: two curves lying in opposite parallel vertical 
planes form the boundary of the surface. We denote the pair of curves either as (zxO(x), zxb(x )) - lying 
in planes perpendicular to the y-axis or as (. zOy(y),zay(y )) if lying in two planes perpendicular to the x- 
axis. 

Let's choose the following base rectangle [(1,2), (5,2,), (5,10), (1,10)]. 


As first example I'll take two parabolas, which should form a surface for 1 < x < 5 and 2 <y< 10. 

2 2 
x - Sx + 11 5 • (x - 6 • x + 5) 

ZXO := , ZXb := - 

2 4 


For representing the two boundary curves in R 3 we have to create space curves. This shall be auto- 
mated: DERIVE shall learn from the appearing variables - x or y - which type is needed. 


C0(kurue, X0 := 0, y0 := 0, Xl := 1, yl := 1) := 
If (UARIABLES(kurue))*1 = x 

SUBSTQx, yO , kurue] , x, x0 + t (x1 - x0)) 

SUBST( [x0 , y, kurue], t, y0 + t(y1 - y0)) 

SUBST([x0, y, kurue], y, y0 + t(y1 - y0)) 

c1(kurue, x0 := 0, y0 := 0, xl := 1, yl := 1) := 
If (UARIABLES(kurue)HI = x 

SUBST( [x , yl , kurue], x, x0 + t(x1 - x0)) 

SUBST([x1, y, kurue], y, y0 + t(y1 - y0)) 

SUBST( [xl , y, kurue], y, y0 + t(y1 - y0)) 


cO ( zxO , 1,2, 5, 10) = [4- 1 + 1 , 2, 8t - 8 t + 3. 
cl (zxb , 1, 2, 5, 10) = [4 • t + 1, IQ, 20 • t • ( 1 - t)] 



Remaining work is not difficult. Pairs of points of the two ’’profile curves” cO and cl, which are gener- 
ated by the same parameter value must be connected. 


fa-fl 

V bj 


•c0 + — -cl oder 


r x \ 

1 -- 

a ) 


cOh cl is the ’’raw form”, assuming one comerpoint of the base 

a 


rectangle being the origin with legths a and b of the rectangle sides. Let's define a general function 
and we don't have to take care for scaling in the future: 


bi_curves ( kO , kl , xO : =0 , yO : =0 , xl : =1 , yl : =1 , zxy ) : =prog ( 
zxy : =IF ( (VARIABLES ( kO ) ) ^l=x , 

( 1- ( y-y0 ) / ( yl-yO ) ) • k0+ ( y-yO ) / ( yl-yO ) . kl , 
(1- (x-xO ) / (xl-xO ) ) • k0+ (x-xO ) / (xl-xO ) • kl , 
(1- (x-xO ) / (xl-xO ) ) • k0+ (x-xO ) / (xl-xO ) • kl ) , 
IF(x0 <x<xlAy0<y< yl , zxy, ?) 



o n 
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Next example shows boundary curves lying in planes parallel to yz-plane. The surface is defined over 
the region described by 0 < x < 6 and 0 <y<8. 


cO(|y - 4|, 0, O, 6, 8), cl (>X( 1 6 - (y - 4) ), 0, 0, 6, 8). 


bi_curues(|y - 4|, 4(16 - (y 


2 

4) ), 0, 0, 6, 8) 



And how does the Tl-calculator perform? 


: b i cur y x ( kO , Kl,xO,xl,yO,yl_} 


: bicurox 
: Func 

: a-<y-yGV<yl_-yG>> + kG+<y-yGV(;yl_-yG>* 
kl 

: EndFunc 


: b i curoy ; kG , k 1 , xQ , y 0 , x 1 , y 1 _ : 

: Func 

: U-Cx-xQVtxl-xG^ + kG+Cx-xGVCxl-xG^k 

: EndFunc 


Up! 1 

■ bicunux] 
x 2 ^9x4 

■ x 2 [Jx4 

|x 2 -Gx + 11 -5-(x 2 - 6-x + 5 ) * 

-^) + 

Done 

... 0-27/2 >-47*v/20+69/4+z5 <x . v >1 

MAIN 

RAD AUTD FUNC i£/30 


Pi^lffi^fcK’cWSi-lPrilolci^ Up! 1 

■ bicunuy(|y - 4| 

JlG -(y - 4) 2 ,G,Q ; 6 ! £ 1 


J-u-(u-B) ly-4il , „ 

x ■ 

6 - 6 J + Iy - 41 

fj-y(y-S) 

ly-4lt , 

6 

6 l + ly 4| + z4(x , y) 


Done 

ans<l>Az4<x,v> 

MAIN RAD AUTD 3D £/30 
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The last generalization is the most difficult one and it finally leads to the Coons-Surfaces. All bound- 
ary lines are - nonlinear - curves. Base figure is again a rectangle and the boundary curves are con- 
nected above the vertices of the base rectangle. 



We produce four parabolas using quadratic regression. 

pari: in xz-plane (y = 0) passing (4,0,1), (0,0,3), (2,0,4) 
par2: in der plane y = 6 passing (4,6,0), (0,6,6), (3,6, -2) 

0 £ 4 


#58 : FIT 


#59 : FIT 


r 2 1 

Lx , a-x + b-x + c J , 


r 2 1 

Lx , a-x + b-x + c_| , 


4 1 

0 3 

2 4 

4 0 

0 6 
L 3 -2 


^ 2 

x 3-x 

= + + 3 

2 2 
AJ 


7-x 37-x 


+ 6 


par3: in plane x=0 passing (0,0,3), (0,2,1), (0,6,6) 
par4: in plane x = 4 passing (4,0,1), (4,3,3), (4,6,0) 
0 ^ y ^ 6 


3-y 7-y 




A 

r 

2 1 

" 0 

3 " 

A 

#60: 

FIT 

Ly i 

a-y + b-y + cj, 

6 

6 






. 2 

1 . 

J 



f 

r 

2 1 

" 0 

1 " 

A 

#61: 

FIT 

Ly i 

a-y + b-y + cj, 

3 

3 






. 6 

0 . 

J 


+ 3 


B-y 3-y 

+ + 1 


IS 


Ans we plot these four curves: 


3-x 


7-x 37-X 


C0 


+ 3, O, 0, 4, 6 


, cl 


+ 6, 0, 0, 4, 6 


, C0 


3y 7-y 


+ 3, 0, 0, 


, cl 


5 • y 3-y 


18 


+ 1 , 0, 0, 4, 6 


4, 6 
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In our first approach we will assume that one vertex of the base rectangle is the origin. Legths of the 
sides of the rectangle are a (in x-direction) and b (in ^-direction). This makes presentation clearer. In 
the original papers only the unitsquare (a = 1,6= 1) is used. 

2 2 2 2 
x 3 ■ x 7 • x 37 • x 3y 7y 5y 

zx0 := + + 3, zxb := + 8, z0y := + 3, zay := + 

2 2 6 6 8 4 18 


3 ' y 

+ i 

2 

[a := 4, b := 6] 

[Z00 := 3, Z0b := 6, zab := 0, za0 := 1] 

The third expression contains the z-values of the ’’corner points”, where the boundary curves connect. 
Later we will automate this, too. 


What follows now is the important and surprising point - for students and for teachers as well. 
The following matrix product delivers the explicit form of the surface: 

y 


X X 

a a 


-z00 -z0b z0y 
-za0 -zab zay 
zx0 zxb 0 


1 - 


2 2 2 
1 6 ■ x • (5 • y - 9) - x • (47 • y + 86y - 432) + 36(3y - 14 y + 24) 


IF 


0<x<a^0<y<b, 


288 


2 2 2 
1 6 ■ x ■ (5 • y - 9) - x ■ (47 • y + 86 y - 432) + 36 (3 y - 14 y + 24) 

288 


So we can proceed plotting the surface: 



Before defining a working function we will try to derive the "Coons-formula". For this purpose we 
need an auxiliary function u(x,y) which decribes the boundary curves. (This part follows the 
references, but we will use DERIVE for the calculations.) 

First of all we add the two different forms of respresentation from bi cu rves () . 

• My H uay . 

a 

It is for sure that this surface does not contain the boundary curves! 

DERIVE enters the stage. 

We take care that ux 0 and uxb are lying on the surface. 


f>- z l 

■ ux 0 + — • uxb + 

( x\ 

1-- 

l b) 

b 

l a) 


u(x, y) 

z(x , y ) := - J ■ u (x , 0) -+- u (x , b) + | 1 - | ■ u (0 , y ) + u (a , y) 

x-u(a, 0) (a - x)-u(0, 0) 

z(x, 0) = u(x, 0) + + 

a a 

As z(x,0) and u(x,0) must coincide we have to subtract the two superfluous expressions from z(x,y) 

f y V x-u(a, 0) (a-x).u(O, 0) 'l y ( * ") * 

z(x, y) := 1 ■ u(x, 0) + u(x, b) + 1 -uCO, y) + u(a, y) 

b J a a Jb LaJ a 

And z(x,b) must coincide with u(x,b), therefore 

x-u(a, b) (a - x)-u(0, b) 


z(x, b) = u(x, b) + 


y ^ f x-u(a, 0) (a - x)-u(0, 0) ^ y f x-u(a, b) (a - x)-u(0, b) 

z(x , y ) := | 1 ■ u (x , 0) + u (x , b) 

b J K. a ajbl. a a 

x 

1-LiCO, y) + u(a, y) 


Let's have a check: 


D 
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#77: [z(x, 0), z(x, b) ] = [u(x, 0), u(x, b) ] 

This does fit!!! 

and what is much more surprising: 

#7B: [z(0, y) , z(a, y)] = [u(0, y) , u(a, y)] 

The other pair of curves are also part of the surface now! 

I rewrite right side of the surface definition using zOO, zxb, zay, zab etc. 

[" y C x-zaO (a - x) ■ zOQ \ y C x-zab (a - x) ■ zOb f x ^ x 

1 - ■ zxO - - + zxb - - + 1 - ■ zOy + zay 

L b J L a a J b L a a J l. a J a 

a-(y-(z00 - zOb - zxQ + zxb) - b- (zOO - zQy - zxO)) - x-(y-(z00 - zOb - zaO + zab) - b- (zOO - zOy - zaO + zay)) 

test :z 

a- b 


and have the question if this expression is identical with the matrix product from above? 


test - 


# 81 : 


1 - 


a a 


" -zOO 

-zOb 

zQy 

-zaO 

-zab 

zay 

_ zxO 

zxb 

0 


y 

1 

b 

y 

b 

1 


= 0 


YES!! 

So we can finish, defining an appropriate function for the general case: 


coons ( zxO , zay , zxb , zOy , xO : =0 , yO : =0 , xl : =1 , yl : =1 , zOO , zOb , zaO , zab , zxy ) : = 
PROG ( 

die z-values in the corner points", 


zOO 

= LIM 

( zxO , 

X 

X 

o 

5 







zOb 

II 

t" 1 

H 

s 

(zOy , 

y, yD 

5 







zaO 

II 

t" 1 

H 

s 

( zxO , 

X, xl) 

, 







zab 

II 

t- 1 

H 

S 

(zay, 

y, yD 

, 







zxy 

= ([1 

- (x 

- xO)/ 

(xl 

- xO) , 

(x - 

xO)/ 

(xl - 

xO) , 

i] • 


[- 

zOO , ■ 

-zOb, z 

Oy; 

- zaO , - 

zab , 

zay; 

zxO , 

zxb , 

o] • 


[1 

- (y 

- yO )/ 

(yl 

- yO) ; 

(y - 

yO) / 

(yl - 

yO) ; 

i]) 

IF (xO < X ; 

< Xl 

a yO < 

y < 

yl , zxy 

, i) 






I'll show some examples to inspire your own experiments. 


2 

sT(l - (x - 1) ), 



2 


|SIN(n x)| 
2 



2 


0 , 0 , 2 , 2 


coons 
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As you can see it works also on the TI-92 / Voyage 200 


fFi-7»10T fit i 

|t Icont-rol 

f'TTV'i 

1X0 

f fht'i 

Uar 

[f ind...l 

[Model 

f 1 

: coonsgerKzxl 

3 ? zay ? zxb ? zGy ? xG 7 yG ? xl ? yl_!> 


: Local zO0 ? zOb , zaO , zab 
: 1 inut-CzxG, x ? xGO+zOG 
: 1 inut-CzGy, y ? yl_)-*zGb 
s limits zxG ,x,xl) *zaG 
: 1 initXzay, y ? yl_)->zab 
: < [ El-Cx-xOWxl-xO), Cx-xGVtxl-xO) , 1 1 ]* 
[ E "Z00, _ zOb ? zOy 1 [ _ za0 7 _ zab ? zay 1 EzxO, zxb 
,0] ]*[h-<y-yOlx<yl_— y0!J ] [ <y-yOV<yi_-y 

&>][i]mi,n 

: EndFunc 


fFiT^OY FT? YTs^Y fht y f? t fs^ t "i 
IfllgebralCalcloiherlFrgnlolciean Up| | 


1 coonsgen 
|sin(n - x)| y 




- J -x(x -2)^1 - -j] ' 


Y'fv 


■ ?vf 1 - —1 + ■ 


“-(V) 


1f~Y 10 


... /2 >+sin <n#y/4>/2-y/4-frsl (x, v>l 

MAIN RAD AUTD 3D E/30 


FtS&H 1 


[^JSlz&^celRe^aphlM^lD^^^l!^ 1 

eyeT=- 4 

xroin=G. 

xroax=2. 

xgrid=lG. 

ynun=G. 

ypiax=2. 

ygrid=lG. 

znun=G. 

znax=l . 

ncontour=G. 






MAIN RAD AUTU 3D 

MAIN RAD AUTD 3D 


A variation of the surface from above: 


coons 


2 

^{1 - (x - 1) ), 



3 |SIN{n x) | 
2 



0 j . 


2 , 


2 
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The base of the next two surfaces is the unit square 0 < x < 1 and 0 < y < 1 . 


#85: coons 


|s-x - 8-x + l| +5 2-7(1 - y) + 1 x -x+1 2-7(1 - y) + 1 


8 


#86: coons 


2 ^ 


2 2 1 2 - y 

2-x -(1 - x) + , 

5 10 


4 1 2 - y 

4 ■ x ■ (x - 1) ■ (2-x - 1) + , 

10 10 



It is not difficult to transfer the surface formula to DP Graph: 



We plot again the surface from above bounded by four parabolas together with a second surface which results of a translation 
of the first one in order to test function coonsQ. 


#87: fll(x, y) coons 


2 2 2 
k 3-x 5-y 3-y 7-x 37-x 

+ + 3, + + 1, 

2 2 18 2 6 6 


3-y 7-y 

+ 6, + 3, 0, 0, 4, 6 

8 4 


#88: f!2(x, y) := coons 


2 2 2 
(x - 1) 3- (x - 1) 5- (y + 3) 3- (y + 3) 7-(x-l) 37-(x-l) 

-+- + 3 , - -+- + 1, - ■+■ 6 , 

2 2 18 2 6 6 


3-(y t J) 7-(y + 3) 


+ 3, 1, -3, 3, 3 
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The Tl-presentation looks a bit poor after the great 
DERIVE and DPGraph plots, but I find it remark- 
able what can be done with these small devices 
algebraically and graphically as well. 



Finally Fd like to check my function once more comparing the function values of the two surfaces 

fll and f 12. 


valuesl := VECTOR (VECTOR (fll(x , y) , y, 0, 6), x, 0, 4) 
values? := VECTOR (VECTOR (f 12 (x, y) , y, -3, 3), x, 1, 3) 
0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 


valuesl - values2 = 


So we don't have to rely only on the optical impression but have an additional numerical confirmation. 


Final comment 

In April I had the occasion to give a 
workshop on Coons-Surfaces. 

The colleagues - teachers from Tyrol 
- enjoyed experimenting and creating 
their own surfaces. So they designed 
their fancy cars, etc. 

My friend Heiner Juen composed two 
surfaces for ’’Josef's Helmet”. Many 
thanks Heiner. I'll keep it in honour. 



References: 

I found the Coons-Surfaces in a Belgian Teachers' Journal [2]. Jan Vermeylen helped providing part 1 [1] of 
this contribution which contains the basics. Many thanks for this together with best regards to Belgium. 

[1] H.van Looy, Coons oppervlakken bij computer aided design , Wiskunde & Onderwijs, 1990 

[2] L. Van den Broeck, Motorkappen en Coons oppervlakken, Wiskunde & Onderwijs, 1996 

This is a related website:: 

http : //olli . inf ormatik . uni -Oldenburg . de/ Graf it i3/ graf iti/ f 1 owl 1 /page 3 . html 

See also additional webpages provided by Ernest Carpenter (User Forum). 
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Introduction to Global Position System 

Two Dimensional Preliminaries 

Carl Leinbach, Bigglerville, USA, leinbach@gettysburg.edu 


NOTE: While working with this lesson, we recommend that you open a 2D-Plot Window and 
then choose the option Window > Tile > Vertically. Later in this lesson, you will be asked to 
modify the range of the window using the Plot Window options Set > Plot Range. Using this 
screen format, you will be able to follow the lesson and view the corresponding graphics displays 
that you will generate. 


Acknowledgement: This activity was inspired by a talk given by Mr. Ralph Irons of Central Shenan- 

3 

doah Valley Regional Governor's School at the 2001 T International Conference in Colombus, OH. 
His presentation can be found in the Proceedings of that Confence call "2001 a Technological Odys- 
sey." The presentation given here is distinctly different from that of Mr. Irons, but his excitement and 
exposition were the inspiration. 


Background - Step 1: 

The classical definition of a circle is: the locus of all points equidistant from a fixed point. The meas- 
ure of the distance from the fixed point is called the radius of the circle and the fixed point itself is 
called the center of the circle. In the cartesian coordinate system a point is designated by a pair of real 
numbers, say (cl, c2). Thus, we need three parameters to completely describe a circle: cl, the x- 
coordinate of the center; c2, the y-coordinate of the center; and r, the radius of the circle. Analytically, 
this tells us that the definition of a circle is the set of all points, (x, y) that satisfy the following equa- 
tion based on the distance formula: 

2 2 2 
#1: (x - cl) + (y - c2) = r 

Derive allows you to plot expressions such as #1 provided the values of the parameters are defined. 

For example, if you wished to see the graph of a circle with center at (-1, 2) and radius 3, you simply 
highlight expression #1, click on the SUB toolbar button and type in the respective new values for cl, 
c2, and r variables (parameters). Click on OK to replace the parameters with your values. Switch to 
the plot window and click on the plot toolbar button or click Insert > Plot. If you receive a message 
that says "Unable to Plot Highlighted Expression", open a 2D-plot window and click on Window > 
Approximate Before Plotting and try again. 

The resulting plot is the set of points, (x, y), satisfying the new values substituted in expression #1. 

What if you know some points on a circle and want to find the equation of the circle? This is the re- 
verse of our previous example. We need to find the values of the parameters cl, c2, and r. Once 
again, expression #1 is our road to success. The first question we need to address is, how many points 
uniquely determine a circle? Since we have three parameters to determine, it is sensible to guess that 
three points should do the job. Of course, the points can not be co-linear. Let's give it a try. 

Suppose we wish to find a circle passing through the three points (1, 1), (-2, 0), and (-1.-1). 


p34 


Carl Leinbach: Introduction to Global Position Ssystem 


D-N-L#58 


First we copy expression #1 into the Derive Command Entry Line by highlighting it and pressing the 
F3 key. Substitute 1 for x and 1 for y. Press the Enter key. The expression will remain in the Com- 
mand Entry Line. Alter the expression to reflect entering -2 for x and 0 for y. Press enter again. Re- 
peat the process for the point (-1, -1). 

This results in the following three equations. 

2 2 2 
#2: (1 - cl) + (1 - c2) = r 


2 2 2 
#3: (-2 - cl) + (0 - c2) = r 

2 2 2 
#4: (-1 - cl) + (-1 - c2) = r 


Derive has a very powerful solver for systems of polynomial equations. 


Press Solve > System. An inquiry box will ask you how many equations you wish to solve. Indicate 
3. A display similar to the one shown below appears. Since you have the equations as part of your 
Derive worksheet, you need only enter the numbers of the expressions that have the three equations as 
shown below. Click on the blank space in the Solution Variables area and then click Solve. 



The following two expressions are displayed on your Derive worksheet. Note that in expression #6 
there are two answers. Why? 


[ 2 2 2 2 2 2 
#5: S0LVE(L(1 - cl) + (1 - c2) = r , (-2 - cl) + (0 - c2) = r , (-1 - 


#6: 


cl) + (-1 - c2) = r J, [cl, c2 , r]) 

1 1 VlO 11 

cl = - a c2 = a r = , cl = - a c2 = a r 

2 2 2 2 2 


yio ' 
2 . 


At this point we enter the solution for cl, c2, and r into expression #1 and plot the circle. 


#7: 



1 l 2 

2 > 


' yio j2 


Exercise 

1 . Choose three arbitrary non colinear points on the plane and determine the equation of a circle 
through these three points. Graph the resulting circle. 
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Background - Step 2 


Now that we have done a brief study of circles, we will consider how they may be applied to GPS 
System. Let's suppose that you are in a flat plane located at the point having grid coordinates (-1, 2), 
but you don't know that fact. Suppose there are transmitters located at grid coordinates, (3, 0), (-3, 1), 
and (0, 1). What you do know is that you are located at distances of the square root of 20, the square 
root of 5, and the square root of 2, respectively from these three points. You will use that fact to 
locate your position on the plane. 


The first thing that you do is write the equations for the three circles centered at the transmitters and 
having radii equal to their distance from your position on the plane. 

2 2 2 
#8: (x - 3) + (y - 0) = 420 

2 2 2 
#9: ( X - -3) + (y - 1) = ^5 

2 2 2 
#10: (x - 0) + (y - 1) = 

Once we have the equations, it is natural to look at their graphs and find the intersection point(s). 
Highlight each of the expressions in the algebra window and plot them. 

Note that even though there are only two independent variables, you need to plot all three of the ex- 
pressions to find your position. Why? 

You can solve the equations exactly using the Derive Solve > System command as we did before. The 
result of this command is shown below. Note that even though we had only two unknowns, we en- 
tered three equations. 

r 2 2 2 2 2 2 2 

#11: SOLVECL (x - 3) + (y - 0) = 420 , (x - -3) + (y - 1) = 45 , (x - 0) 

2 21 

+ (y - 1) = n J, [X, y]) 


#12: [x = -1 a y = 2] 

At this point, let's review what you know about circles. In the following exercises draw a graph to 
illustrate that your answer is correct. 

Exercises: 

2. Is it possible to find cl, c2, and r so that expression #1 has exactly one point, (x, y) that will 
satisfy the equation. 

3. Is it possible for two distinct circles to intersect at no point? one point? 
two points? 

4. Is it possible for three distinct circles to intersect at no point? one point? two points? 

5. Why is it not possible to find three or more distinct circles that intersect in exactly three points? 
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A "2 Dimensional” GPS 


NOTE: The velocity given for the signal sent from the transmitters has no relation to the actual speed 
with which signals are sent from GPS sattelites or radio transmitters. We are using this velocity only 
to keep the magnitude of the numbers we are dealing with more manageable so that the numbers do 
not get in the way of the exposition. In part II we will deal with a more realistic situation. 

Prior to beginning this unit, resize your graphics screen so that the range on the x-axis is (-300, 900) 
with 4 intervals, and the range on the y-axis is (-200 , 800) with 5 intervals. Also, because accuracy is 
important in GPS calculations set the decimal accuracy to 10 digits. Use the Algebra Window's De- 
clare > Simplification Settings to set the Precision Digits to 10. 

You know that you are located at a point somewhere in this plane region. Assume that scale on the 
graphics screen is in miles. You activate your GPS to locate your position. It receives signals from 
three transmitters that have been placed in the plane at locations of locl=(75, 150), loc2=(510, 110), 
and loc3=( 150,20). The transmitters all know their own position and the exact time as determined by 
an atomic clock in the transmitter. The transmitter clocks are synchronized on a regular basis so that 
the time at each transmitter is absolutely accurate. Each transmitter broadcasts its coordinates and the 
time. The broadcast signal travels at 200 miles per second. 

Your GPS is able to receive the signals and record the time that each signal is received. However, 
your GPS does not have a clock that is synchronized with the transmitters. Thus, there is a time time 
discrepency, say d, between the time recorded by your GPS receiver and the the actual time as deter- 
mined by the synchronized atomic clocks in the transmitters. The GPS records the signals that it re- 
ceived. The times for each transmitter are subtracted from the time your GPS received them. These 
are recorded as tl, t2, and t3. The numbers correspond to the times from the respective transmitter 
locations. . 

#13: tl := -0.8139056099 

#14: t2 := 0.7980845608 

#15: t3 := -0.4685532484 

The actual times that it took to travel from the trasmitters to your GPS are tl+d, t2+d, and t3+d, where 

d is the time discrepency of your GPS. Thus, we have three unknown values: the x and y coordinates 
of your position and d, the time discrepency of your GPS. 

Your distance from each of the transmitters is the velocity of the signal multiplied by the time that it 
took to travel to the GPS. 

#16: dl := 200 -(tl + d) 

#17: d2 := 200-(t2 + d) 

#18: d3 := 200-(t3 + d) 
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But these distances are the radii of circles centered at each transmitter and passing through your posi- 
tion on the plane. We are now in the situation of Step 2 above. The one difference is that we now 
have three unknowns: x, y, and d. However we have three equations to solve. It may work out. 




2 

2 

2 

#19: 

(X - 

75) 

+ (y - 150) = 

dl 



2 

2 

2 

#20: 

(X - 

510) 

+ (y - 110) : 

= d2 



2 

2 

2 

#21: 

(X - 

150) 

+ (y - 20) = 

d3 


We simply use the Solve > System option as shown above and then approximate the result to obtain 
the following result. NOTE: the first result is not shown here. It is the exact symbollic solution to the 
system. It is rather messy since it involves expressions with sums of radicals and fractions of VERY 
large numbers. We want numbers that we can recognize. However, this will be a decimal approxima- 
tion and we may loose some accuracy. 

r 2 2 2 2 2 2 

#22: SOLVE(L(x - 75) + (y - 150) = dl , (x - 510) + (y - 110) = d2 , (x - 

2 2 21 
150) + (y - 20) = d3 J, [x, y, d]) 

#23: [x = 452.5761121 a y = 96.9671036 a d = -1.092505996, x = 132 a y = 145 

a d = 1.1] 


Because we have three unknowns, we, once again, have two possible solutions to the system of equa- 
tions. In this case we can rather easily eliminate one of them. Note that if d = -1.092505996 then each 
of the distances, dl, d2, and d3 will be negative. This is nonsense. We would have received the sig- 
nals before they were sent. Thus, you can conclude that your position is at (132, 145) and the time 
discrepency is 1.1 second. Define d := 1.1 and graph expressions #19, #20, and #21. 

NOTE: This can be done by holding down the <Ctrl> key and highlighting all three expressions fol- 
lowed by switching to the 2D-plot window and clicking on the plot button. 

#24: define d:=l.l by editing this line 
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After viewing the graph clear the value for d by typing d := . This action frees the symbol d so that it 
can be a variable. 

#25: clear d by editing this line to read d := 

It is time to try one on your own. 

First clear the graphics screen by choosing the Edit > Delete All Plots from the 2D-plot menu. 

Assume that you are working with a different GPS receiver, so you have no idea of the time difference 
between the transmitter clocks and your GPS clock. You are still in the same region of the plane as in 
the example. Here we go! 

Exercise 

6. Your GPS receives signals from three transmitters located at positions in the plane given as (320, 
190), (50, 100), and (570, 50), respectively. The time differences recorded by your GPS receivers 
and the three transmitters are given as tl := 3.113941029 seconds, t2 := 3.113941029 seconds, and t3 
:= 4.102775637 seconds, respectively. Enter these times by editing the lines below. 


#26: 

enter 

ti 

by 

editing 

this 

1 i ne 

#27: 

enter 

t2 

by 

edi ti ng 

thi s 

1 i ne 

#28: 

enter 

t3 

by 

edi ti ng 

thi s 

1 i ne 


Edit the next three lines to give the formulas for the distance from each of the three transmitters. Don't 
forget that you need to include the time discrepency between the transmitters and your GPS in these 
formulas. 


#29: 

enter 

the 

di stance 

to 

transmi tter 

#1 

by 

edi ti ng 

thi s 

1 i ne 

#30: 

enter 

the 

di stance 

to 

transmitter 

#2 

by 

edi ti ng 

this 

1 i ne 

#31: 

enter 

the 

di stance 

to 

transmitter 

#3 

by 

edi ti ng 

this 

1 i ne 


Now that you have the distances, you can enter the three equations that you will need to solve in order 
to determine your position.on the plane. Edit the three expressions below so that they contain the 


equations. 

#32: 

equation 

#33: 

equation 

#34: 

equation 


After editing the above lines, solve the equations you entered using the Solve>Equations option. In 
this case, the distances for either of the answers are positive. However, we can eliminate one of the 
aswers rather easily. It lies outside of the region of the plane where we know we are located. Write 
the value for d, the time discrepancy, and plot the three circles. 

#35: value for d here 

As a final exercise, you will demonstrate your understanding of the process by constructing a problem 
on your own for others to solve. 
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Exercise 

7. Choose a set of 4 points that lie within the 2-D plot region of the plane that you have laid out and 
write down their coordinates. One of these points will be your position on the plane. The other three 
will be the locations of the three transmitters. Determine the distance of your point from each of the 
transmitters. Finally decide on a value for the variable, d, which is the time discrepency between the 
transmitters clocks and the clock on your GPS. 

a. Determine the time that it takes for each of the signals sent from the 
transmitters to reach your GPS. 

b. Using d, determine the values for tl, t2, and t3, the differences 
between the time sent by each transmitter and the time the signal 
is received. 

c. Write out your problem in a manner similar to the example and 
problem given above. 

d. Solve the problem and compare the answer with the coordinates 
for your position. This is a check that you correctly constructed 
the problem. 

#36: InputMode := Word 


One more letter from Ernest Carpenter: 

Ernest Carpenter (South Africa) 

The article on Turtle Graphics on the TI-92, in DNL57, also brought back some nice memories when I 
first obtained my TI-92+ (before that I was actually a die hard HP48GX fan© and also used the QLIB 
referred to in the “Actuarial Math” article) and TG was one of the first programs I load onto it! 

Now that I’m finished babbling, just a few quick questions to get your option on a matter or two: 

> Do you think we will see a scaled down version of Cabri3D for the TI-89/92+/V200? 

> Do you think from an education point of view whether the Tl Graphing Calculators (or Graph- 
ing calculators in general TI,HP, Casio) would still be around for a while or will it disappear 
and make place for more windows based handhelds like the PocketPC or PalmOS devices? 
The main reason for this question is to see if developers or users in general should still focus 
on developing software/apps for the current devices or shift their efforts towards the new 
technology marked? 

> I know that the use of technology in RSA, especially on a secondary (high school) mathemati- 
cal level, is almost non-existent. How do the students/pupils in Europe benefit from using 
technology, like for instance the TI-92 or Derive? 

Finally, any indication when you will ‘publish’ the articles “Tools for 3D-Problems” and “CAD-Design 
with DERIVE and the Tl”? 

(I know I’m persisted with the above mentioned question, but maybe my persistence would pay-off©) 
Thanks again for all the information and articles, I appreciate it! 

Enjoy your day! 

Regards 

Ernest Carpenter 

CADD Co-ordinator 
Strategic Projects 
De Beers(NM) 
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Titbits from Algebra and Number Theory (30) 

by Johann Wiesenbauer, Vienna 

Let's consider the following innocent looking problem: Given a positive integer n, 
is there a right triangle with rational sides a,b,c and area n? In other words, we 
are looking for rational solutions (a,b,c) of the system consisting of the two 
equations: 

a 2 +b 2 =c 2 , ^ = n (*) 

Since (a,b,c) is a rational solution w.r.t. to the area n, if and only if (ta,tb,tc) is a 
rational solution w.r.t. to the area t 2 n for any nonzero integer t, there is no loss 
of generality if we always assume that n is squarefree in the following. A 
squarefree positive integer n, for which the system of equations above has ra- 
tional solutions, is called a congruent number. 

Just to get a feeling for the problem, let's look now into some small cases. For 
example, are there rational solutions for n = 5 and n = 6, respectively? Well, this 
shouldn't be too hard for n = 6, as the "Egyptian triangle" with the sides 3,4,5 
clearly solves the problem. On the other hand, it might take you a while to see 
that n = 5 is a congruent number and e.g. a = 3/2, b = 20/3, c = 41/6 is a rational so- 
lution of (*) here. Just in case you wonder why I didn't start with the first 
three squarefree natural numbers, namely n = 1, 2, 3, the answer is very simple: 
There aren't any rational solutions in those cases, but this is far from being ob- 
vious. (E.g., the case n= 1 was settled by Fermat around 1650.) 

There is a completely different way of looking at this problem, which will give us 
valuable insights though. Using Derive the following computations show that for 
a fixed n we can assign to each rational solution (a,b,c) of (*) a rational point (x,y) 
with y^O on the curve y 2 = x 3 -n 2 x and also the other way round. 

2 2 2 2 2 
a- (a - -/(a + b )) a ■ (-/(a + b ) - a) a-b 

x := , y := , n := 

2 2 2 . 

2 3 2 

y -x + n -x = 0 


[x :=, n :=] 







D 


-N-L#58 


Johann Wiesenbauer: Titbits 30 


p41 


Algebraic curves of the form y 2 =f(x) , where f(x) is a polynomial over some field 
F without multiple roots, are called elliptic curves. In many applications, it suf- 
fices to consider elliptic curves of the special form 

y 2 =x 3 +ax + b (a,beF) (**) 

where the condition that the polynomial on the right-hand side has no multiple 
roots becomes the condition that its so-called discriminant A = 4a 3 + 27b 2 doesn't 
vanish in F. 


In the special case F=[R we can take advantage of sliderbars for the parameters 
a and b (e.g. in the range [-10,10]) to produce the following pictures: 
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Note that only the first two curves are elliptic curves, while the other two are 
not due to singularities. 

Maybe you have already heard of elliptic curves and their usefulness, e.g. in solv- 
ing Diophantine equations (the most notable example is A. Wiles' proof of FLT!) 
and in encrypting messages. Hence, unless you are already an expert in this field, 
you might wonder: What is so special about them? 

Well, there are many possible answers, but one of the most outstanding fea- 
tures of elliptic curves is certainly the fact that it is possible to define an addi- 
tion on their points which turns them into abelian groups. 

Considering again the field F = R, the basic procedure to obtain the sum of two 
points P and Q is the following: Draw the chord through P and Q, find its inter- 
section point R with the curve and define P+Q as the mirror image of R w.r.t. the 
x-axis, as shown in the picture below. 


The rough description above must be supplemented though in the case, where P 
and Q have the same x-coordinate. Then either P and Q coincide, in which case 
the chord becomes a tangent or/and P and Q are mirror images to each other 
w.r.t. the x-axis. In the latter case, we use the usual convention from projective 
geometry that the intersection point R with the curve (and hence also P+Q) is 
the point O at infinity. In particular, this point O has to be added to the points 
of the elliptic curve and plays the important role of the neutral element for our 
addition. 

Let's now use the informal description above to compute the coordinates of P+Q, 
where P = (x 1 ,y 1 ) and Q = (x 2 ,y 2 ), under the assumption that P + Q + 0. Then 


P 
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is the slope of the chord or tangent, respectively, and after an easy computation 
we get for the coordinates of P + Q = (x 3 ,y 3 ) the equations 

x 3=k 2 - x i- x 2> y 3 =-Yi +k( x i - x 3> 

It is important to notice that these formulas, actually derived under the condi- 
tion F = R, also make sense for an arbitrary field F, and although there is no 
longer a geometrical interpretation then, it is still true that the set of points of 
an elliptic curve along with this addition is always an abelian group. 

When trying to implement this addition (and, of course, this is the next point on 
our agenda!) we must somehow specify the field F. In order to avoid making 
things to complicated for a first approach (for example, by allowing arbitrary 
fields for which the basic operations are given by functions defined by the 
user), I will consider in the following only the case, where F is either a subfield 
of the field C of complex numbers or a residue class ring mod p for some prime 
p. Then a routine for the addition of two points U and V on an elliptic curve given 
by an equation (**) could look like this: 

"inv(x, p) := 

If p = 0 
1 /x 

INVERSE_MOD(x , p) 

add(u, v, a, p := 0 , k_) := 

Prog 

If L 4 I + V 4 .I = CO 

RETURN [1/Cl/uil + 1M1) , 1/CL/L42 + 1M2)] 

If L 4 I = vj,l 

Prog 

If M0D(m2 + m2, p) = 0 

RETURN [co, co] 

k_ := M0 D(O , uj, 1 a 2 + a) ■ nnv(2 ■ uj,2 , p) , p) 

If k_ = ? 

RETURN GCD(2.i42, p) 

Prog 

k_ := M0D((vj,2 - uj,2) ■ inv(vj,l - ml, p) , p) 

If k_ = ? 

RETURN GCDCv^l - ml, p) 
v := M0D(k_ A 2 - ml - ml, p) 

[v, M0D(- m2 + k_'(ml - v) , p)] 


A few comments might be in order here. First, you should know that in the case 
p = 0, i.e. if F is assumed to be a subfield of C, those reductions modp will do no 
harm, although they are of course superfluous! Second, in the case F = Z p , we 
will occasionally p allow to be a composite number, which leads to an "elliptic 
pseudocurve", for which addition is sometimes undefined. In those cases, a non- 
trivial divisor is returned as the result of the addition, which could be just 
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the thing you were after. 

For the time being, we're dealing with a more down-to-earth problem though, 
namely the determination of some points on an elliptic curve again given by an 
equation y 2 =x 3 +ax + b. If F = Z p , and p is not too large, we should be able to 
compute even all its points. As for other fields, we restrict ourselves here to 
the most important case F = Q, where the following routine computes all points 
up to a certain user-defined "height" h, which is an upper bound for |r|,|s| in a 

representation x = ^ of the x-coordinate of a point (x,y) on the curve. At any 

rate, the point O at infinity is added. 

points(a, b, p := 0, h := 100, s := +1, h := 0 , f_, u_ := {[<», «>]}, x_ := 0, y_) := 

If p > 0 
Loop 

y_ := SQUARE_R00 T(x_ a B + a-x_ + b, p) 

If NUMBER? (y_) 

u_ :z u_ u {[x_, y_], [x_, M0D(-y_, p)]} 

x : + 1 

If x_ = p 
RETURN u_ 

Loop 
h_ :+ 1 
WRITE(h_) 

If h_ > h 
RETURN u 

f_ := SELECT (GCD(n_, h_) = 1, n_, 0, h_)/h_ 

Loop 

x := FIRST (f_) 

y_ := ^(x_ a B + a-x_ + b) 

If RATIONAL? (y_) 

u_ := ADJOIN ([x_, s-y_], u_) 
y_ := VC- x_ a B - a-x_ + b) 

If RATIONAL? (y_) 

u_ := ADJOIN ([-x_, s-y_], u_) 
x_ := l/x_ 

y_ := 7 (x_ a B + a-x_ + b) 

If RATIONAL? (y_) 

u_ := ADJOIN ([x_, s-y_], u_) 
y_ := 7(- x_ a B - a-x_ + b) 

If RATIONAL? (y_) 

u_ := ADJOIN CC-x_, s-y_], u_) 
f_ := REST (f_) 

Iff_= [] ex it 

Before applying this routine, one should make sure though that the parameters 
a and b are chosen appropriately, i.e. A = 4a 3 + 27b 2 must not vanish in F, as men- 
tioned above. This can be accomplished by 

B 2 

regular?(a, b, p := 0) := M0D(4-a + 27 -b , p) t 0 


At last, we are ready to apply the tools we have gathered so far to our 
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introductory problem! (Note that all the curves y 2 =x 3 -n 2 x at issue are non- 
singular due to the nonzero discriminant A = 4a 3 +27b 2 =-4n 6 in F = Q.) Okay, 
maybe still not quite, as it will be convenient to have a routine that translates a 
rational point U*0 on those curves with y^Ointo the corresponding right tri- 
angle with rational sides a,b,c. 


u - n 2-n-u u 

1 11 


tnangle(u, n) := 


po-ints(-36 , 0) 

" I" 6 720 ' 

-[-6, 0], [-3, ±9], [-2, ±S], , ± , [0, 0], [6, 0], 

v L 49 343 . 

TABLE(tnangle(u , 6), u, SELECT(u t 0 a u t ®, u, points(-36, 0))) 


25 35 " 

, ± , [12, ±36], [IS, 

4 S . 


±72] , [™ , m ] \ 


[-3, ±9] 

[-2, ±S] 

6 7 


[3, 4, 5] 

[4, 3, 5] 

120 7 1201 


49 343 J L 7 10 70 

25 35 1 f 7 120 1201 


4 S J L 10 7 70 


[12, ±36] 
[IS, ±72] 


[3, 4, 5] 
[4, 3, 5] 


As you can see above for n = 6, by searching for rational points up to the default 
height h=100 / we have got another (far less obvious!) right triangle with rational 
sides a = 7/10, b = 120/7, c = 1201/70 and area 6. 

What about the other squarefree numbers n in the range l<n<10? Take a look 
at the following table with all the rational points found up to height 100: 


TABLE(po-ints(-n „ 0) , n , SELECT (SQUAREFREE(n) , n „ 1 , 10)) 
Dnsp]ayFormat:jCompressed 

r i {[-i 


{[- 1 . 0 ] , [ 0 , 0 ] , [ 1 , 0 ] , [®,®]} 
{[- 2 , 0 ], [ 0 , 0 ], [ 2 , 0 ], [»,«]} 
{[-3,0], [0,0], [3,0], [®,®]} 
100 1 f 25 


9 27 

6 720 

49 343 



" 25 

75 " 

. [0, 

0], [5,0], 

,+ 


. 4 

S . 

,[0, 

" 25 

0] , [6,0], 

35 ' 


. 4 

S . 


LS,±72], [®,®]l 


' r 63 735 1 r 49 1176 1 

[-7,0], ,± , ,± , [0,0], [7,0], [25, ±120], [®,«] 

, L 16 64 J L 25 125 J 


{[- 10 , 0 ], [ 0 , 0 ], [ 10 , 0 ], [®,®]} 
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Well, I already told you that 1,2,3 aren't congruent numbers and from the table 
we can conclude that not only 5 and 6, but also 7 is a congruent number for sure. 
But what about n = 10? It is true that we haven't found rational solutions of 
y 2 = x 3 -lOOx except for the trivial ones, which are of no use here, but the tan- 
talizing question remains: Was this because our search limit was simply too low 
or can it be really proven that n = 10 is not a congruent number? 

If this is what you asked yourself (and hopefully you did!) then I have good news 
and bad news for you. First the good news: There is a remarkably simple test by 
J.B.Tunnell (1983), in order to find out whether a squarefree positive integer is 
congruent or not. The bad news: Only "no" answers are provably reliable, while 
"yes" answers are not. It must be said though that no counterexample has been 
found so far, in fact, such a counterexample would contradict the famous con- 
jecture by Birch and Swinnerton-Dyer (one of the 7 millennium problems, each 
worth $ 1,000,000 for the first solver), which is widely accepted by the commu- 
nity of mathematicians. 

For this test, we must compute the numbers r and s of integer solutions of the 
equations x 2 + 2ay 2 + 8z 2 = — and x 2 + 2ay 2 + 32z 2 = — , where a = 2-(nmod2), and 

ci ci 

simply check whether r = 2s or not. The following routine will do the trick. (By 
the way, if you feel like a programming challenge, I cordially invite you to have a 
try at this task yourself before looking at my solution below. Although it is not 
too hard in principle, one has to take care not to get lost in too many cases and 
subcases, resulting in ugly if-statements.) 

congruent?(n , a_, u_, r := 0 , s_ 0, x_, y_, u_, v_) 

Prog 

:z 2 - M0D(n, 2) 
n : / a_ 

u :z REST (SORT ({0 FLOOR (» } ■ {0 FLOOR (./(n/ (2 ■ a_))) })) 

Loop 

If U_ = [] 

RETURN r_ = 2-s_ 
x_ :z FIRST (FIRST (0) 
y_ :z FIRST (REST (FIRST (u_))) 
v_ :z x_ a 2 + 2 ■ a_-y_ A 2 
z_ -/((n - v_)/8) 

If INTEGER?(z_) 

Prog 

r_ :+ (2 - 0 A x_) ■ (2 - 0 A y_) ■ (2 - 0 A z_) 

If EVEN?(z_) 

s_ :+ (2 - 0 A x_) ■ (2 - 0 A y_) ■ (2 - 0 A z_) 
u := REST (u_) 

As you can see from the table below, 10 is not contained among the numbers be- 
low 100, that passed the test, which is a proof that it isn't congruent, indeed. 
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TABLE(congruent?(n) , n 

, SELECT (SQUAREFREE(n) , n 

, l, 

100))’ 





' 1 

2 

3 

3 

6 

7 

10 

11 

13 

14 

13 

17 

19 

21 

22 

23 

26 

29 

. false 

false 

false 

true 

true 

true 

false 

false 

true 

true 

true 

false 

false 

true 

true 

true 

false 

true 

30 

31 

33 

34 

33 

37 

38 

39 

41 

42 

43 

46 

47 

31 

33 

53 

57 

58 

true 

true 

false 

true 

false 

true 

true true 

true 

false 

false 

true 

true 

false 

true 

true 

false 

false 

39 

61 

62 

63 

66 

67 

69 

70 

71 

73 

74 

77 

78 

79 

82 

S3 

S3 

86 

false 

true 

true 

true 

false 

false 

true 

true 

true 

false 

false 

true 

true 

true 

false 

false 

true 

true 

87 

89 

91 

93 

94 

93 

97 












true 

false 

false 

true 

true 

true 

false _ 













From a closer look at the table you will also learn that all squarefree numbers of 
the form 8k+5, 8k+6 or 8k+7 in the given range passed the test and hence are 
congruent numbers, as it has been proven by others that all results of Tunnell's 
criterion are correct, if n < 1000. In fact, congruent numbers not of this form 
are a minority, like 34, 41, 65 in our example above. 

Before coming to an end, let's also talk a little bit about the case, where F = Z p 
for some prime p, which plays an important role in modern cryptography. One of 
the most important issues there is the computation of the number N P of points 
on an elliptic curve E again given by the equation y 2 =x 3 +ax + b (a,beF). It turns 
out that the algorithms used for this computation are very similar to those in- 
troduced in my last column to solve the discrete logarithm problem (DLP) and it 
will soon become clear why. 

A very first, though very crude implementation, would make use of our routine 
points (a, b,p) and simply count the number of resulting points, i.e. 

Np(a, b, p) DIM(points(a, b, p)) 

As the following simple example shows this is far too slow though. 

regular?(l, 1, 1009) = true 

NpCl, 1, 1009) = 1034 (18.4s) 

A better idea is to compute for every x in the range {0,l,...,p-l} not the points 
with this x-coordinate, but only the number of possible y-values, i.e. the number 
of solutions of y 2 = x 3 +ax + b modp. According to Euler's criterion from the the- 
ory of quadratic residues this number is given by l + ((x 3 +ax + b) (p_1)/2 modp), 
which yields the implementation 

3 (p - l)/2 

Np(a, b, p) := p + 1 + X(M0DS(0< + a-x + b) 

NpCl, 1, 1009) = 1034 (0.03s!!) 


p), x, 0, p - 1, 1) 
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On the other hand, the following computation 

6 

Np(l, 1, NEXT_PRIME(10 )) = 1000727 

still took 156.5s on my PC! It looks as if we are not yet arrived at a method for 
counting points on elliptic curves modp, where p has say 10-20 digits! 

Hence, let's use some results from the theory of elliptic curves, which may bring 
us further. In the first place, there is the well-known theorem by Hasse, which 
states that the number of points of an elliptic curve E mod p, which we denote by 
#E P in the following, always lies in a relatively close vicinity of p + E To be more 
precise, the absolute difference |#E p -(p + l)| can never get larger than 2 A /p . 

How can we exploit this feature for our purposes? 

Well, as you might know, the order of an element of a finite group is always a 
divisor of the group order. This means that starting with an arbitrary point U on 
the elliptic curve E p , there is always some multiple kU with k "nearby" p + 1, such 
that kU = O. If the order of Ef is "big enough", more precisely not smaller than 
the length 4^/p of the "Hasse interval" [p + l-2 A /p,p + l + 2 ^/p], then this k must 
be #E P itself, as it is the only candidate. 

Okay, let's do the programming now. In the first place, we need a routine that 
computes for any integer n the multiples n\J of a point U on the elliptic curve E p , 
which is given as usually by an equation y 2 =x 3 +ax + b. Again, try to write this 
routine yourself, before looking at my solution below. 

mult(u, n, a, p, s_ := [«>, «>]) := 

If n < 0 

mult([u 4 ,l, - 142 ], -n, a, p) 

Loop 

If n = 0 
RETURN s_ 

If 0DD?(n) 

s_ := add(s_, u, a, p) 
u := add(u, u, a, p) 
n := FLOORCn, 2) 

We will later have many opportunities to see the "power" of this small jewel of 
programming, hence I'll skip examples here. 

The next thing we must do is to solve the equation kEJ = O for a given point U, 
where the integer variable is the unknown here. This is where our algorithms for 
the DLP come into play. (Hopefully, you can see this yourself, although the op- 
eration is written additively, whereas I used the multiplicative notation to intro- 
duce the DLP in my last column!) For example, by slightly adapting the 
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baby-step giant-step algorithm we get the following routine. (Note that the 
variables u, a, p have the same meaning as above!) 

B5G5(u, a, p, j_, i_ := 0, r_, s_, t_, v_, x_) := 

Prog 

r_ := CEILING(p A (1/4)) 

s_ := mult(u, p + 1 - 2-r_ A 2, a, p) 

t_ := mult(u, 2 ■ r_ , a, p) 

v_ := ITERATES (add (u , w_, a, p) , w_, [«>, <»], r_) 

x_ := v_ COL 1 

Loop 

If MEMBER?(FIRST(s_), x_) 

Prog 

j_ := POSITION(FIRST(s_), x_) - 1 
If s_ = v_j,(j_ + 1) 

J- := -j- 

RETURN p + 1 + 2.(i_ - r_) ■ r_ + j_ 
s_ := add(s_, t_, a, p) 

i_ :+ 1 

Finally, to get the order of U, we must successively remove from the returned k 

all primes q for which still — U = 0 holds. 

q 

ord(u, a, p, o_, q_) := 

Prog 

o_ := BSGS(u, a, p) 

q_ := (FACTORS (o_)) COL 1 

Loop 

If q_ = [] 

RETURN o_ 

Loop 

If mult(u, o_/FIRST(q_) , a, p) t [°°, °°] exit 
o_ :/ FIRST (q_) 

If M0D(o_, FIRST (q_)) > 0 exit 
q_ := REST(q_) 


Well, are you ready for the most advanced Np(a,b,p) so far? Here you are! 

Np(a, b, p, g_ := 1 , h_, s_ := 1, x_, y_) := 

Prog 

If p < 229 

RETURN p + 1 + I(M0DS((x A 3 + a-x + b) A ((p - l)/2), p) , x, 0, p - 1, 1) 
Loop 
g_ :+ 1 

If JACOBI (g_, p) = -1 exit 
Loop 
Loop 

:= RANDOM (p) 

s_ := JACOBI (k_ a 3 + a-x_ + b, p) 

If s_ — ±1 exit 
h_ := g_ A ( (1 - s_)/2) 

y_ := SQUARE_R00T(h_.(x_ A 3 + a-x_ + b) , p) 
h_ ord([M0D(h_-x_, p) , H0D(h_-y_, p)], M0D(h_ A 2-a, p) , p) 

If h_ > 4-4p 

RETURN (p +!)■(! - s_) + s_-CEILING(p + 1 - 2- Jp, h_)-h_ 





p50 


Johann Wiesenbauer: Titbits 30 


D-N-L#58 


For those of you, who, who wonder how this works in detail, here comes a little 
support. At first, you should be able to see that s_ = 1 means that for a ran- 
domely chosen x in {0,l,..,p-l} there exists a point (x,y) on the curve indeed. In 
the opposite case s_ = -1, there is still no need to despair, as every curve E p has a 
"twin brother" E’ p , called its quadratic twist, with the property that there is a 
point with a given x-coordinate on E p if and only if there is no such point on E’ p . 
Its equation is given by 

y 2 = x 3 +ag 2 x + bg 3 

where g is some quadratic nonresidue modp. From this easily follows that 

#E P + #E’ P = 2(p + 1). 

Hence, if we know one of the two summands on the left, we know also the other 
one! This will do for all primes p>229 as noted by Mestre. If this condition is not 
fulfilled, then we simply use one of the previous formulas, which are fast enough 
for those small primes! 


Okay, I have to come to an end, as Josef is already tearing out his hair, both 
because of the length of this paper and its late delivery. Hence, let me conclude 
with some small examples for the routines above. First, let's take a really huge 
number k, say the lcm(l,2,.. ,30000), to form the multiple kU for U = (1,1) mod F 7 . 
(After all, Fermat told us that F 7 is a prime, didn't he?) 

DIM(LCM([1 30000])) = 13013 


r 

multL[l, 1], LCM([1, 


7 


2 

, 30000]), 139, 2 + 


1 J 


- 59649589127497217 122.7s!!) 


Whew, what's that? Obviously Fermat was mistaken and our routine has re- 
ported a 17 -digit prime factor p of his number. And now the acid test for our 
third try at Np(a,b,p), which is the time-consuming part of ord(u,a,p)! 
ord([l, 1], 139, 59649589127497217) = 29824794689700405 1106.7s) 


FACTOR (29324794689700405) = 3 ■ 5 ■ 37 -419 -1021 -4783 ■ 26263 


It worked! And what's more, it worked in a reasonable time! What you saw above 
is the basic principle of the celebrated ECM (=Elliptic Curve Method), intro- 
duced by H.Lenstra Jr. in 1987 . As you see, by a fortunate coincidence the or- 
der of the point ( 1 , 1 ) on the elliptic curve E29824794689700405 splits up into many 
small prime factors, the largest being 26263, which is smaller than our bound 
30000 above. Actually, without this prime factor, the order would be even 
"smooth" w.r.t. to the much lower bound 4783 (or say 5000). In an actual 
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implementation of the ECM, one takes advantage of this fact by introducing two 
bounds Bi and B 2( where Bi covers all prime factors with one possible exception 
whereas B 2 (usually 10 to 100 times larger than Bi) covers also this one-off. 

As for the following Derive-implementation, here is the meaning of the vari- 
ables: n is the number to factor, U is a point on the elliptic curve used, a is the 
coeff icient in the equation of the curve, and s and t are the two bounds above. 
(Note that, as many times before, again the seemingly missing coeff icient b of 
the curve equation is determined implicitly by the coordinates of the point Ul) 


ECM(n, u, a, s, t, !_:=[], m_ := 0, u_, p_ := 2, q_) 
Prog 
Loop 

u := mult(u, p_ A FL00R(LN(s) , LN(p_)), a, n) 

If NUMBER?(u) 

RETURN u 

p_ := NEXT_PRIME(p_) 

If p_ A 2 > s exit 
Loop 

u := mult(u, p_, a, n) 

If NUMBER?(u) 

RETURN u 

p_ := NEXT_PRIME(p_) 

If p_ > s exit 
u_ := mu"lt(u, p_, a, n) 
q_ := NEXT_PRIME(p_) 

Loop 

If NUMBER?(u_) 

RETURN u_ 

If p_ > t 
RETURN 1 
Loop 

If q_ - p_ < m_ exit 
m_ :+ 2 

v_ := mult(u, m_, a, n) 

If NUMBER (v_) 

RETURN v_ 

1_ := APPEND(1_, [v_]) 
u_ := add(u_, l_4,((q_ - p_)/2), a, n) 

P- := q_ 

q_ := NEXT_PRIME(q_) 


Will this be notably faster than our crude computation above? Sure! 


r 7 

2 

ECMl2 + 1 


[1, 1], 139, 5000, 30000j 


= 596495 891 2 749 7 21 7 (4.89s!) 


Don't take this time as the total running time for ECM though, as there are 
quite a few failures until you find the "appropriate" value a = 139 above. Okay, 
okay, let's stop here ... 


(As always, questions or remarks to i.wiesenbauer@tuwien.ac.at ) 
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Recommended Readings 

Liebe Kolleginnen und Kollegen 

Hansruedi Vollmer und ich haben beschlossen, die 2. ueberarbeitete und verbesserte Aufla- 
ge unseres Buches "Skalarprodukte, Schwingungen, Signale" auf dem Bildungsserver 

www.swisseduc.ch 

kostenlos anzubieten. Bitte weitersagen. 

Mit den besten Gruessen 

Hansruedi Vollmer und Hansruedi Schneebeli 

Hansruedi Vollemer and Hansruedi Schneebeli offer a free download of their book “Scalar 
Products, Oscillations, Signals” (in German) from www.swisseduc.ch . 


Information provided by Guido Herweyers, Belgium 

You can find the pdf-version of “A CASE FOR CAS” online at 

http://www.uhasselt.be/scholennetwerk/iq/CAS/CAS.html 
(from www.t3vlaanderen.be publicaties ) 

This is a T 3 -Europe product written by Ian Forbes, Gert Schomacker, Guido Herweyers, 
Rene Hugelshofer and Josef Bohm. 


Information provided by Robert McCollum, USA 

This is just a reminder to send in a proposal to speak at the NCTM Chicago Regional to be 
held in Chicago Sept. 20-22, 2006 at the Hyatt McCormick Place. The deadline for submissi- 
ons is JULY 1 1 , 2005 (this summer!) 

I would love to see you sharing your expertise with CAS at this conference. 

The URL is www.nctm.org/speak/ 

They are particularly looking for CAS speakers but most anything you submit would be ap- 
preciated. 


Take Care, 
Bob 


